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WALL INTERFERENCE IN A TWO-DIMENSIONAL-FLOW WIND TUNNEL, WITH 
CONSIDERATION OF THE EFFECT OF COMPRESSIBILITY 
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SUMMARY 

Theoretical tunneh-waU corrections are derived jor an airfoU 
of finite thickness and camber in a two-dimensional-flow wind 
tunnel. The theory takes account of the effects of the wake of 
the airfoil and of the compressibility of the fluid, and is based 
upon the assumption that the chord of the airfoil is small in 
comparison with the height of the tunnel. Consideration is 
given to the phenomenon of choking at high speeds and its 
relation to the tunntmrnU corrections. The theoreoretical 
resuUs are compared with the small amount, of low-speed ex~ 
perimental data available and the agreement is seen to be satis- 
factory, even for relatively large values of the chord^height ratio. 

INTEODUCnON 

The need for reliable wind-tunnel data for the design of 
high-performance aircraft has led in recent years to attempts 
to make the conditions of the tunnel tests conform more 
closely with the conditions prevailiag in flight, especially with 
regard to the Eeynolds and Mach numbers. Because of 
practical limitations in. size and power, most existing wind 
tunnels, whether high speed or low speed, are not capable 
of providing full-scale Reynolds nmnbers for all flight con- 
ditions. In order to obtain the highest Reynolds numbers 
possible imder the circumstances, it is necessary to use 
models dimensions of which are as lai^e as possible relative 
to the cross-sectional dimensions of the tunnel test section. 
The effect of such large size is to make the test conditions 
depart further from the conditions prevailing in flight by 
increasing the magnitude of the tunnel-wall interference. In 
the case of tests at high Mach nmnbers, the interference is 
increased still further by the tendency of the flow pattern 
of a compressible fluid, if unrestrained, to expand as the 
Mach mimber of the imdisturbed stream increases. Since 
the walls of a closed-throat timnel restrain certain of the 
streamlines at a flxed distance from the model, this expansion 
is prevented, and the tunnd-waH interference and corrections 
become progressively larger as the Mach number increases. 
The results obtained in the tunnel must therefore be corrected 
accmrately for the effects of wall interference if they are to 
be appKed with confidence to the prediction of free-flight 
characteristics. 

In tests at high Mach numbers an additional complication 
arises. The effect of a model ia a closed-throat tunnel may, 
in a sense, be thought of as equivalent to that of a constric- 
tion ia the throat of the tunnel. The resulting convergiog- 



diverging nozzle formed by the model and the tunnel walls 
then has roughly the same characteristics at high speeds as 
the usuoj supersonic nozzle; that is, for some Mach number 
less than unity ia the undisturbed stream, sonic velocity is 
reached at all poiats across a section of the tunnel at the 
position of the model, and the flow ia the diverging region 
downstream of this section becomes supersonic. When this 
occurs, increased power input to the tunnel has no effect 
upon the velocity of the stream ahead of the model, the 
additional power serving merely to increase the extent of the 
supersonic region ui the vicinity of the model. At this 
point the tunnel is said to be "choked" and no fm-ther 
iacrease in the test Mach nimiber can be obtained. The 
value of the Mach nxunber at which choking occms is thus of 
extreme importance, siace it determines the upper limit of 
the range of Mach nmnbers which can be obtained with a 
given combination of model and tunnel. 

In testing airfoils to obtaiu section characteristics at 
subsonic speeds, it has become common practice in modem 
closed-throat wind tunnels to have the model span the tiumel 
80 that supporting struts and their accompanyiag iater- 
ference effects are entirely eliminated. If the tunnel has a 
cross section of rectangular shape, this arrangement results 
in a flow which is essentially two-dimensional. 

The wall interference for such a two-dimensional-flow 
wind timnel has been the subject of nmnerous investigations, 
the results ia general beiag expressible as series in ascendiag 
powers of (c/h), where c is the chord of the airfoil and h the 
height of the tunnel. The effect of wall interference upon 
the flow of an ideal fluid about a symmetrical airfoil at zero 
angle of attack is determined to the order (c/hy by Lock ki 
reference 1 and by Glauert ia reference 2. The interference 
for an infinitesimaUy thin, cambered airfoil at a small angle 
of attack in an ideal flxiid is given by Glauert to the order 
(c/hy ia reference 2, and iavestigations for the special case 
of a flat plate have been carried out to a higher order of 
accuracy by several writers. While the present report was 
being prepared, work by Goldsteia appeared (references 3 
and 4) ia which the interference is determined to the order 
(cfli)* for a general cambered airfoil of finite thickness ia an 
incompressible fluid, no restriction being made in the general 
results as to the magnitude of the camber, thickness, and the 
force coefficients. A still later paper by Goldsteia and 
Young (reference 6) gives the modifications necessary in the 
previous results to allow for the effect of fluid compressibility 
to the order (c/hy. 

1S5 



156 



REPORT NO. 782 — NATIONAL ADVISORT COMMITTEE FOR AERONAUTICS 



In the present paper, the tunnel-waH corrections are 
determined to the order (c/A)* for the general airfoil in 
a compressible fliiid for Mach numbers below that at which 
chokiiig occurs. It is assumed that the thickness and 
camber of the airfoil are small and that the interference 
velocities are everywhere small as compared with the velocity 
of the undistiu'bed stream. A discussion is also included of 
the Mach number at which choMng occurs. The various 
results presented are of essentially the same nature as those 
which already have appeared separately in the references 
cited, but the methods of development and certain of the 
final results axe different, especially with regard to the 
interference associated with lift. The validity of the final 
corrections is examined by comparison with the available 
experimental data. The equations also are compared with 
the results of references 3, 4, and 5, and the afore-mentioned 
differences are discussed. 

The discussion is limited to airfoils placed midway between 
the upper and lower walls of the tunnel. Mathematical 
symbols are defined as introduced in the text. For reference, 
a li^ of the more important symbols and their definitions is 
given in appendix B. 

DEVELOPMENT OP CORBECTION EQUATIONS 

In an analysis of tunnel-wall interference it is desirable 
to look upon the theoretical development of the tunnel-waU 
corrections as consisting of two parts. First, it is necessary 
to determine the manner and extent to which the tunnel walls 
alter the field of flow about the airfoil from what it would 
be if they were not present. Second, it is necessary to 
calculate the effect of these alterations upon the measured 
characteristics of the airfoil. The development of the 
correction equations of this report has been divided into 
these t^vo general sections. 

In reference 6, the use of the method of superposition to 
determine the pressure distribution over the surface of an 
airfoil section in free air is presented. It is shown that in 
the calculation of the flow at the surface of a thin airfoil of 
small camber, the effects of camber and thickness may be 
considered independently. This follows directiy from the 
fact that the velocities induced by the vortex sheet used to 
represent camber and those induced by the source-sink system 
used to represent thickness are simply additive in their 
effect on the flow over the airfoil. 

To treat the problem of wall interference, it is again con- 
venient to consider the thickness and camber effects sep- 
arately. The flow changes associated with airfoil thickness 
are foimd by considering the interaction between the tunnel 
walls and the base profile of the airfoil, the base profile 
being defined as the profile the airfoil would have if the 
camber were removed and the resulting symmetrical airfoil 
placed at zero angle of attack. The interference effects 
associated with airfoil camber are found by analyzing the 
interaction between the timnel walls and an ind&nitesimally 
thin airfoil having the same camber as the actual airfoil. In 
addition to the interference effects associated with airfoil 
thickness and camber, it is necessary to consider a fmUier 
alteration of the field of flow caused by the confining influence 
of the timnel walls upon the airfoil wake. When the indi- 



vidual effects promoted by the interference betewen the walls 
and the airfoil thickness, camber, and wake are known, the 
total alteration in the flow at the au'foil is foimd by super- 
position, and the characteristics of the airfoil in the altered 
field of flow are compared with the characteristics in free air, 
This comparison leads to skaple formulas which enable the 
prediction of the free-flight characteristics when the charac- 
teristics in the tunnel are known. 

The method of superposition, which is fundamental to the 
entire analysis, is in general inapplicable to compressible flow 
as the differential equation for such flow is nonlinear in the 
physical plane. The separate solutions which are super- 
posed are obtained, however, by assuming that the ai rfoil is 
of small thickness and camber and that the induced velocities 
are thus small as compared with the velocity of the undis- 
turbed flow. On the basis of this assumption the equation 
of compressible flow becomes a linear differential equation — 
namely, Laplace's equation (references 7, 8, and 9) — so that 
superposition of velocities is, in this case, technically per- 
missible. Furthermore, the tunnel-wall corrections are in 
most cases rather small relative to the experimental quanti- 
ties being corrected, so that it is not thought that the use of 
this approximate method will lead to large errors in the final 
corrected quantities. 

INFLUENCE OF TUNNEL WALLS UPON FIELD OF FLOW AT AIRFOIL 

Thickness effect. — The interaction between the base pro- 
file and the waUs of a two-dimensional-flow tunnel has been 
considered by Lock for the case of an incompressible fluid 
(reference 1 ; a discussion of Lock's method is also given by 
Glauert in reference 2). Lock's method of analysis is 
essentially to introduce an infinite series of images of the 
base profile such as to satisfy the condition that there is no 
flow normal to the walls, to replace each image by a suitable 
source-sink doublet, and to calculate the velocity induced 
at the base profile by this system of doublets. It is shown 
that the net effect of the tunnel waUs upon the flow at the 
base profile is to increase the effective axial velocity of an 
incompressible stream by the amount 

(AjFOi=A<rF' (1) 

where 

V apparent stream velocity at airfoil as determined from 

measurements taken at a point far ahead of model 
ff a factor dependent upon size of airfoil relative to tunnel 
A a factor dependent upon shape of base profile 
The factor <r is defined by the equation 

where (c/h) is the ratio of the airfoil chord to the tunnel 
height. The factor A can be determined for any base profile 
jrom the relation 

where 

yi ordinate of base profile at chordwise station x 
dyjdx slope of siirface of base profile at a; 
Pf^ base-profile pressure coefficient at x in an incom- 
pressible fluid 
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(Ifc will be noted that the quantity \P in references 1 and 2 is 

equivalent to ^ Ac? in the notation of this report.) Yalues of 

A for a number of base profiles are given in table I. 

In appendix A, it is shown for linear theory that the effect 
of compressibility upon the streamwise induced velocity 
at a given point a large distance above or below a body ia a 
uniform stream is such as to multiply the velocity increment 
for incompressible flow by the factor AfJ"^ where M 
is the Mach number of the flow far upstream from the body. 
Applying this result to the velocity induced at the base 
profile by each of the airfoils in Lock's system of images, it 
can be seen that for a compressible fluid the increase in the 
effective axial velocity in the tunnel is 

AiF^= [l_(^.)^j3/» A^y- (4) 

where M' is the apparent Mach number; that is, the Mach 
number corresponding to the velocity V. The result of 
equation (4) has also been obtained by an independent 
procedure in reference 5. 

The value of A, as given by equation (3), is actually 
slightly larger than that which would be obtained using the 
first-order theory of this report. Moreover, second-order 
theory would probably indicate that the value of A in the 
compressible case would be slightly modified for the effect 
of compressibility. A corresponding term for the drag 
correction due to streamwise pressure gradient has been 
shown in reference (10) to be slightly modified by the influ- 
ence of compressibility. It is dearly inconsistent then to 
employ the second order value of A without including the 
effect of compressibility upon its value. This additional 
complication, however, has not been introduced since the 
velocity correction due to thickness is of major importance 
in those cases for which (fijh) is large and it is customary, and 
properly so, to employ lai^e values of (c/A) in low-speed 
investigations only. 

Consideration of the symmetry of the base profile and of 
the system of images used by Lock to simulate the effects 
of the tuimel walls indicates that the interaction between 
the walls and the base profile does not induce velocities 
normal to the center line of the tmmel. Similarly the base 
profile does not affect the longitudinal velocity gradient in 
the tunnel at the position of the airfoil. 

Wake effect. — ^In the wake of an airfoil moving through a 
real fluid, the total head of the fluid is less than in the region 
outside the wake. This reduction arises from the increase 
in thermal energy caused by fluid friction m the boundary 
layer and in the wake itself and by any shock waves which 
may exist in the vicinity of the airfoil. Considering a section 
normal to the wake, it may be said that the static pressure 
across the stream is nearly constant if the section taken is 
not too close to the trailing edge of the airfofl. It follows 
that the reduction in total head which exists within the wake 
must appear almost entirely as a decrease in the local dynamic 
pressure of the fluid. This decrease arises primarily from 
a reduction in the local velocity and secondarily from the 
reduction in local density which accompanies the increased 
temperature within the wake. Thus, since the local velocity 

8-13110 — BO ^11 



and density within the wake are both less than in the external 
flow, the mass-flow rate per unit area is less inside the wake 
than outside. This condition prevails both in the tuimel 
and in free air. In the tunnel, however, the requirement of 
continuity of flow between a transverse section upstream 
from the airfoil and a section across the wake necessitates, 
in addition, that the mass-flow rate per xmit area outside of 
the wake is greater than the mass-flow rate per miit area 
ahead of the airfofl. In order to satisfy this requirement, 
the velocity in the tunnel outside of the wake must be 
greater than that of the undistinrbed stream. This fact 
implies that as the flow proceeds down the tuimel the velocity 
of the main portion of the stream imdei^oes a gradual 
increase from the value prevailing in the undisturbed 
stream ahead of the model to some higher value downstream 
of the airfofl. This does not hold true in free air, where the 
velocities of the main flow upstream and downstream of the 
model are equal. The interference between the wake and 
the tunnel walls thus gives rise at the position of the model 
to a velocity increment and a velocity gradient which are 
not present in an unlimited stream. Further, as required 
by BemouUi's equation, the velocity gradient is accompanied 
by a longitudinal pressure gradient which likewise would not 
exist in free air. 

To determine the magnitude of these effects the procedure 
is briefly as follows: Two stations in the tunnel are considered, 
one far upstream from the model and one far enough down- 
stream so that the wake has spread to the walls and the 
velocity is again uniform across the tunnel. The difference in 
static pressure between these two stations is evaluated as a 
function of the measured drag of the airfoil. The pressure 
gradient at the airfofl can be related to this pressure differ- 
ence and hence to the drag of the airfofl by a convenient 
analytic device, which is essentially the some as that used 
by Groldstein (reference 3). The airfofl and its wake are 
considered to be replaced by a fluid source located at the 
position of the airfofl. It is specified that conditions far 
upstream in the resulting hypothetical flow must be the same 
as those existing in the actual stream. With this provision, 
the magnitude of the velocity and static pressure far down- 
stream can be determined as fimctions of the upstream 
conditions and the strength of the soxirce. The strength is 
then related to the drag of the airfofl by requiring that the 
static pressure difference promoted between the two stations 
in the tunnel by the source flow is the same as that which 
actuaUy exists when the airfofl and wake are present. The 
tunnel walls can then be replaced by an infinite system of 
such sources directly above and below the position of the 
airfofl at intervals equal to the height of the txmnel. The 
system of image sources alone, however, would induce a 
small finite negative velocity at infinity upstream, so that 
it is necessary to superpose on the flow field an additional 
uniform flow of equal velocity in the positive direction in 
order to satisfy the original requirement that the conditions 
far upstream shafl be unchanged. The velocity of this 
flow, which is readfly determined as a function of the source 
strength and hence of the airfofl drag, then gives the velocity 
increment caused at the airfofl by the interference between 
the wake and the waUs. The longitudinal velocity and 
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pressure gradients at the position of the airfoil are found in 
terms of the drag by evaluating the flow induced at that 
point by the image sources. It is apparent that this entire 
method of analysis fails to satisfy the actual condition as 
regards the velocity at infinity downstream. This dis- 
crepancy arises out of the fundamental difference between 
the actual flow in the wake and the source flow by which it 
is represented and is unavoidable as long as this representa- 
tion is used. 
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Fiarms 1.— Soorce s^ntem for anal^ of wake effect. 

Consider the flow in a closed two-dimensional-flow wind 
tunnel, as shown in figure 1. At a station far upstream, 
the effect of the model upon the flow is negligible, so that 
the velocity V, the density p', the static pressiu-e p', and 
the absolute temperature T' are constant across the stream. 
At a station far downstream, where the wake has spread to 
the walls, the velocity V", the density p", the pressure p", 
and the absolute temperature T" are again constant across 
the stream. 

The difference between the pressures p' and p" can be 
related to the measured drag of the airfoil by means of the 
conditions of continuity, conservation of energy, and im- 
pxilse and momentiun, together with the state relations for 
a perfect gas. The condition of continuity is given by 



(5) 



and, if it is assumed that the flow is an adiabatic process, 
conservation of total enei^ requires that 



or 



2gJc,T' (T'' \\\ 



(6) 



where. 

g gravitational acceleration 



J mechanical equivalent of heat 

Cp specific heat of gas at constant pressure 

In modem wind tunnels the walls of the test section are 

flared slightly to compensate for the growth of the bomdary 

layer on the walls, and only the drag of the airfoil therefore 

need be considered. The impulse-momentum equation can 

be written 



h 



=p'-p''+p'iy'y-p"{V"Y 



or 



w~ p'iV'Y p' \V') 



(7) 



where D' is the drag of the airfoil and d' the drag coefficient 
referred to the apparent dynamic pressxure g'. 

The velocity of sound Ve in the undisturbed stream is 
related to the absolute temperature by 



(F/)*=7i23"=(T-l)fl'Jcpr' 



(8) 



where y is the ratio of specific heats and B is the gas constant. 
By means of this relation, equation (6) con be written 



and, from equation (5), 



(?-0 



(t-i)(MO' 



(9) 



£__ 



(7-l)(M0= 



(10) 



The state equation for a perfect gas then provides the relations 

-i 



p" _ p"T' _T" 
p'~p'T'~T' 



(y-1)(M')' 



(11) 



and 



P' _ RT _ (.Vc'y _ 1 



(12) 



Substitutions from equations (9), (10), (11), and (12) into 
equation (7) gives 



rpir 

1 — — 

J- IJl/ 



(T-1)(M')' 



+1- 



(7-l)(M')' 
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from which it can be found that 



I (Y-l)(M')j (y+lHMr 



^[l-(MOT-T(MO»^(02+>y(Mo{2-^(f)]} 



(13) 



For airfoils usually employed, the factor is smaJl. 

Expanding the above expression and neglecting terms con- 
taining c/ to powers higher than the first gives 



^i+T'TWrrcJS^ ^^^^ 



By means of this relation, together with equations (9) and 
(10), the static pressure difference in terms of g' is obtained 
from equation (7) as 

a' — \h)v+i-(M'y\ ^^^^ 

Now, consider the airfoil and wake to be removed and 

replaced by a fluid source of strength Q. If the flow con- 
ditions far upstream are maintained unchanged, the mass 
flow far downstream is then 



or 



(16) 



where the subscript « denotes conditions now prevailing at 
the latter station. For reversible adiabatic flow 



i(MO.[(^)-i]j 



Since it is to be expected that iV,"/V') will be dose to unity, 
the right-hand side of this equation may be expanded in 

ascending powers of j^^-^^ — ij, and terms containing 

powers higher than the first neglected. Thus, 



^=x-<fl'[(¥)'-] 

and equation (16) becomes 



(17) 



~rni- s— 



It is reasonable to assimie that the ratio Ji^^^ is small as 



compared with unity. The solution of the preceding equa- 



tion to the first order in 



Q 
h/V' 



IS 



(18) 



Bernoulli's equation for reversible adiabatic flow can be 
written 

^=7t^[-(-^<-).[(^)-x])*] 

Since ( -yrj is close to unity, this may be replaced by the 
approximate relation 

Substitution from equation (18) and neglect of the term 
involving the square of then gives 



2' ~hyV l-{M'y 



(19) 



Comparison of equations (15) and (19) shows that the source 
strength required to promote the same pressiure difference as 
actually arises from the confinement of the airfoil wake is 



Q 



1+ 



(20) 



The tunnel walls are now replaced by an infinite system 

of sources of strength Q spaced h distance apart and located 
directly above and below the position of the airfoil as shown 
in figure 1. This image system together with the souirce 
which has been placed at the position of the airfoil satisfies 
the requirement that the flow at the plane of the tunnel wall 
shall be tangential to the wall. 

As shown in the first of equations (All) of appendix A, a 
source of strength <? in a uniform flow of compressible fluid 
will induce at a distance r from itself a streamwise velocity 



AV-- 



- <2 f 
2irpr I ^l—AP (l—AtP an* 



cos <i> 



*)1 



where <f> is the polar angle of the point in question and p 
and M are the density and Mach number of the undisturbed 
stream. By virtue of this relation, the streamwise velocity 
A^y induced at a point on the center line of the tunnel by 
the entire system of image sources is 



Q 



m-1 vp r„ 



cos ^, 



Vl-(W [1 - (M'r mi' <l,J 



where r„ and ^„ are the radial distance and the polar angle 
of the point relative to the source a distance mh above or 
below the center line and p' and M' are the density and 
Mach number of the undisturbed flow in the tunnel. If the 
distance from the position of the airfoil to the point on the 
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center line is denoted by x (taken positive downstream), this 
equation can be "written 



A»V'= 



Q 



or 



■irX 



(21) 

It can be seen by setting z= — oo in equation (21) tbat the 
image sources induce at an infinite distance upstream a 
velocity 

(AaT^O— = -2p'h[l-(M')'] 

In order to satisfy the original requirement that conditions 
far upstream remain, unchanged, this velocity must be 
counterbalanced by the superposition of a imiform flow of 
equal magnitude but opposite sign. The addition of this 
flow at all points in the field will result in a speeding up of 
the general flow at the position of the airfoil by the amoimt 



Q 



2p'h[l-{M')'\ 
or, substituting the source strength from equation (20), 

If the factor r is defined as 



(22) 



the velocity increment induced at the position of the air- 
foil by the interference between the wake and the walls may 
thus finally be written for air (7=1.4) as 



(23) 



The longitudinal velocity gradient produced at the position 
of the airfoil by the flow from the image sources can be foimd 
by differentiating equation (21) with respect to x and then 
setting a;=0. This gives ficnaUy 

dV'_ d(A^V') _ -kQ 
dx~ dx ~WW^^WW^ 

or, by virtue of equation (20), 

dV fl+(7-l)(M0^ 



dx 



■' fl+ 
"=("[1 



■ (MOT"" 



7rc/V'e 



(24) 



It already has been noted that the interference associated 
with the Uiickness of the airfoil has no effect upon the 
longitudinal velocity gradient at the position of the airfoil. 
It will be seen later that this also is true of the interference 
associated with airfoil camber. Equation (24) thus gives 
the total velocity gradient for the complete airfoil and wake. 



The total pressure gradient at the position of the model 
then is given by Bernoulli's equation as 



dy 
ax 



— -ST 



or, substituting from equation (24) and setting 7=1.4, 



dv_ \ l+QA{M'Y \ TrCa'(i'c 
l[i-(M')T''l"6F~ 



(25) 



It is apparent from the symmetry of the system of image 
sources that at the center line of the tunnel the interference 
between the wake and the walls has no effect upon the ve- 
locity normal to the direction of the stream. 

It is shown later in this report that the camber of the air- 
foil does not affect the stream velocity at the airfoU. Equa- 
tions (4) and (23) together thus give the total increase 
m. velocity for the complete airfoil and wake. The effective 
or true velocity V at the model is therefore 

^=4+ [i-(i^0T- ^-+ 'wMl-^' -^''1 

It is evident that a correction to the apparent velocity in a 
compressible flow implies corrections also to the apparent 
density, dynamic pressure, Reynolds number, and Mach 
number. These corrections are readily obtained on the basis 
of the usual assumption that the flow is adiabatic. It is as- 
sumed that the correction terms ore small as compared with 
unity, so that squares and products of these terms may bo 
neglected. 

The true density p at the model is connected with the 
apparent density p' by the isentropic relation 

P=,'(i-^(Mt[(^)-i]P (27) 

Substitution from equation (26) gives, after expansion as an 
ascending power series and neglect of correction teions higher 
than the first order. 



P=P' 1 



(MO' 



[i-(M')T'' 



Ao- 



{M'y[l+QA{M')^ , 
I- {My 



(28) 



The true dynamic pressiure 2=2''^ ^ related to the 
apparent dynamic pressm:e g' by the equation 



By means of equations (26) and (28) this can be written to 
the first order as 



2=2 y-^[i—{M'yf'* 



[2-(M0^ [l-i-0 4(M0'] 

\-{M'y '^"^ 



(29) 



The true Reynolds number R is given in terms of the 
apparent Reynolds UTunber R' by the equation 
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where n and n' are the coefficients of viscosity coixesponding 
to V and V. According to von Kdrmdn and Tsien (refer- 
ence 11), tlie coefficients of viscosity are related to the cor- 
responding absolute temperatures by 

y.'\T') 

For reversible adiabatic flow it can be shown that 

T=T' [l -"^{M'Y [(t^J- l]) (30) 

which after substitution from equation (26) becomes to the 
first order for air (7=1.4) 

OA(M'y 0.4(M0'[1+0.4(M01 ,) . 

By means of these relations together with equations (26) 
and (28), the true Reynolds number may be written 

n--R> U , 1-0;7(M0» . . [1-0.7(3/01 [1-H0.4(M0^ ^ . ,1 

(32) 

The true Mach number M is related to the apparent 
Mach number M' by the equation 



where Ve and VJ are the velocities of sound corresponding to 
p^and v. Since the velocity of sound in a gas is directly 
proportional to the square root of the absolute temperature 
alone, this equation may also be written 



M=M'(^,)(f 



With the aid of equations (26) and (31) the true Mach nmn- 
ber then may be written to the first order 



M=M' 



^^[l-(M')1^/»^'^^ 



[1 +0.2(3/01 [1-1- 0.4(^/01 ^, 
1-(M0' ' 



') 



(33) 



At low Mach numbers, the terms containing to/ in the 
correction equations are usually negligible as compared with 
the terms containing Ao-. At supercritical Mach numbers, 
however, where the drag coefficient is very large, the terms 
with rCtf' are predominant. 

Numerical values of the compressibility factors appearing 
in equations (26), (29), (32), and (33) are given in table 11. 

Camber effect. — The theory* of the infinitesimally thin, 
cambered airfoil in free air is developed by Glauert in refer- 
once 12 (pp. 87-93). In this development the camber line 
is replaced by a sheet of continuously distributed, bound 
vortices. The flow induced at any point on the camber line 
by this system of vortices is obtained by integration and is 
combined vectoriaUy with the flow of the undisturbed stream 



to give the direction of the resultant flow. The distribution 
of vorticity is then determined from the condition that the 
resultant flow at all points on the camber line must be tan- 
genital to the camber line. 

In the actual calculation of the induced velocity, it is 
assumed that the vortices may be distributed along the 
chord line rather than along the camber line and that the 
induced velocity at any chordwise station on the camber 
line is the same as the induced velocity on the chord line at 
the same station. If the origin of coordinates is taken at 
the leading edge of the airfoil (fig. 2), with the positive 
X-axis along the chord line and the positive i/-axi8 directed 
upward, the induced velocity {v)i in an incompressible fluid 
at any point Xo on the chord line is 




FlouBE 2.— Mean-camber line In bee air. 



■ dx 



(34) 



where dVldx is the vorticity per imit length at the point 
X and c is the chord of the airfoil. The direction of this 
velocity is normal to the x-axis. 

Glauert (reference 7) has shown that a first approxima- 
tion to the velocity induced at any point by a simple vortex 
in a compressible stream can be obtained by simply multi- 
plying the velocity induced at the same]^oint in an incom- 
pressible stream by the factor 



l—M* sin» <^ 



(35) 



where Mis the Mach number of the undisturbed flow and «^ 
the polar angle of the point in question as measured from 
the direction of flow of the imdisturbed stream. For points 
on the chord of an airfoil which is inclined at a small angle 
to the direction of the undisturbed stream the polar angle 
^ is small, and the factor (35) is sensibly equal to 

If it is assumed that the effect of a vortex sheet in a com- 
pressible fluid may be obtained by superposing the effects 
of elementary vortices, the velocity induced at any point 
Xo on the chord line in a compressible fluid is 



v= 



dT J 



2fl- 



Jo X— Xo 



(36) 



If the undisturbed velocity of the free stream is taken 
equal to the true velocity V at the airfoil in the tunnel, 
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the condition that the resultant flow shall be tangential to 
the camber line requires that, for all points on the airfoil, 



= j — a 



'7~dx 



(37) 



where dyjdx is the slope of the camber line at a^o and a is 
the true angle of attack; that is, the angle the xmdisturbed 
stream makes with the diord line in free air. (See fig. 2.) 



(a) 




FioDBB 3.— Mean-camber line In tnnnd. 



The problem of the infinitesunally thin, cambered airfoil 

in a two-dimensional-flow tunnel can be investigated by the 
method of images; that is, the effect of the upper and lower 
walls of the tunnel can be simiilated by introducing an 
infinite lattice of alternately inverted but otherwise identical 
image airfoils above and below the original airfoil, as shown 
in figure 3(a). By this artifice the direction of flow at the 
position of the upper and lower walls can be made to coincide 
with the plane of the walls, which is the required condition 
of flow. As in Glauert's analysis of the airfoil in free air, the 
camber line of the airfoil and of each of its images is replaced 
by a sheet of continuously distributed vortices, the vortex 
distribution of all sheets being identical in magnitude but 
alternately reversed in sign. The flow induced at any point 
on the camber line of the original airfoil by the entire vortex 
system is then obtained by integration. As before, the dis- 
tribution of vorticity must be determined so that the result- 
ant of the induced velocity and the stream velocity is 
tangential to the camber line of the airfoil. 

For the detailed calciilation, the coordinate system is 
taken as shown m. figme 3(b). The origin of coordinates is 
taken on the center line of the tunnel at the leading edge 
of the airfoil. The positive as-axis extends downstream par- 
allel to tihe undisturbed flow, and the positive y-axis is 



directed upwards. It is assumed that the vortices may be 
distributed along the a;-axis and the induced velocities calcu- 
lated at points on this axis. This arrangement is somewhat 
different from the system employed for the airfoil in free 
air, where the x-axis was taken along the chord line; how- 
ever, since the angle of attack is assumed to be small, the 
difference is of no consequence. 

It is evident from figm-e 3(b) that, for an airfoil midway 
between the upper and lower walls of the tunnel, the axial 
velocity induced at any point on the x-axis by any one ' 
image is nullified by the velocity induced by the corres- 
ponding image on the opposite side of the timnel. It 
follows that airfoil camber does not affect either the true 
axial velocity or the longitudinal pressmre gradient in the 
tunnel at the position of the model. 

The vertical velocities induced at any point on the ai-axis 
by any one image and its counterpart are, however, addi- 
tive. Thus, for corresponding images situated at mh and 
—mh, respectively, the vertical velocity {i/„^i induced at 
the point in an incompressible fluid is 



(»'rJ*=2(-l)"JJ- 



dV . ( ^ A 
W sm [^^r.—^) 



dx 



2irr, 



or 



dV 

1 C'~dx (*~"^»)^ 

{Vrm)i=- (-1)"J„ (a;_ai,)*+(mA)» 



(38) 



where dV'Idx is the vorticity per unit length at the point x 
in the timnel. 

It win now be assimied that the chord of the airfoil is 
reasonably small in comparison with the height of the tunnel. 
This being the case, the approximation 

is sufficiently precise for purposes of this analysis, and 
the term (i— a%)* in the denominator of equation (38) may 
be neglected. 

The vertical velocity (!)',)< induced by all the images is 
then found by superposition as 



(-1) 



-Jo ^ 



dx 



12K- 



dx 



(39) 



This equation can be corrected for the effect of compressi- 
bility by means of expression (35). If, as was assumed, the 
chord of the airfoil is reasonably small as compared Avith the 
tunnel height, the polar angle (j) of any point Xg on the airfoil 
with respect to any point x- on an image is nearly a right 
angle, so that in this case the factor (35) is sensibly equal to 



The vertical velocity induced in a compressible stream by 
aU the images is then 
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The vertical velocity v't, induced at a point cKb by the 
vortex sheet belonging to the airfoil itself, is given by 
equation (36) if F' and v't are substituted for r and t>, 
respectively. 

The total vertical induced velocity »' at any point on 
the airfoil in the tunnel is then the sum of v'r and v't; that is, 

The condition that the resultant of the induced velocity 
and the true axial velocity at the airfoil shall be tangential 
to the camber line requires tiiat, at all points on the camber 
line, 



v'_dyc , 



(42) 



where a' is the angle of attack of the airfoil in the tunnel; 
that is, the angle the chord line makes with the center line 
of the tunnel. The true velocity y rather than the apparent 
velocity V is used in equation (42), since the vortex system 
used to represent the cambered airfoil in the tunnel is actually 
operating in a stream of velocity "Fwhen the airfoil thickness 
and wake are present. 

BELATIONS BETWEEN CHABACTEBISTICS OF AIRFOIL IN TUNNEL AND 

IN FBBE AIR 

The preceding sections provide the basic information 
required for tiie development of relations between the charac- 
teristics of the airfoil in the tunnel and in free air. The 
relations for the lift and moment coefficients and angle of 
attack are derived from the equations of the preceding 
section by on extension of the method of Fourier series 
employed in Glauert's theory of thin airfoils (reference 12, 
pp. 87-93). To this end, the vorticity distributions for the 
airfoil in the tunnel and in free air are each represented by a 
trigonometric series, the two series being similar in form but 
having undetermined coefficients. By means of the equa- 
tions of the preceding section, general relations are found 
between the coefficients of the two series. These general 
relations ore then specialized to meet the requirement that 
the airfoil shall have the same value of the cotangent term of 
the series in the tunnel and in free air, this requirement being 
shown to be necessary to assure that the essential character- 
istics of the pressure distribution will be sensibly the same 
in both cases. By means of the relations between the 
coefficients, expressions are then derived for the lift and 
moment coefficients and angle of attack of the airfoil in free 
air in terms of the characteristics measured in the tunnel. 
The corresponding drag coefficient in free air can be found 
from the drag measured in the tunnel by subtracting the 
pressure drag caused by the interference between the walls 
and the wake and referring the remaining drag to the true 
instead of the apparent dynamic pressure. Knally, a 
method is presented for correcting airfoil pressure distribu- 
tions for the effect of tunnel walls. 



To carry out the analyBis, points on the airfoU are defined 
by a new coordinate 9 such that 



and 



a;=2c(l— cos 6) 



dx=2 emid dd 



(43) 



(44) 



The distribution of vorticity along the chord of the airfoil 
in free air is represented, after Glauert, by the trigonometric 
series 

^=2V^AoCot^d+f^nsmne^ (45) 

Equation (36) then gives the induced velocity at any point 
8 on the airfoil as 



— -4o+S -^n COS tW 
n-l 



and equation (37) for the slope of the mean-camber line 
becomes 



J*=a- ■^f^^Ao+ ^/1=WJ^^A, COS nO 
The coefficients are then given by the relations 



(46) 



^cosnBde 
;o ox 



(47) 



For the airfoU in free air the coefficients A„ for n^l are 
thus functions of the camber-line shape only and are inde- 
pendent of the angle of attack. The coefficient is a 
ftmction of both the camber-line shape and the angle of 
attack. 

The chordwise lift distribution in free air is given by 
dL „dr 2 dr 

which after substitution from equation (45) can be written 
in coefficient form as 



(48) 



Equation (48) illustrates the well-known fact that in free 
air the chordwise lift distribution consists essentially of 
two distinct parts. The one part, contributed by the sine 
terms and generally referred to as the basic lift distribution 
(reference 13), depends in magnitude and form only upon 
the shape of the mean-camber line. The other part, defined 
by the cotangent term and referred to as the additional lift 
distribution, is fixed in form and depends in magnitude upon 
the angle of attack as well as upon the camber-line shape. 
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The distribution of vorticity for tJie airfoil in the tunnel is 
represented by 



rlV f 1 " 



dx 



(49) 



Substitution of this expression, together with expressions 
(43) and (44), into equation (41) gives, after int^ation, 

[^.'-jZ^ (2^'+^i')] COS <?+ 

SA.' COSTlfll 
n=2 J 

where o- is as defined by equation (2). Equation (42) for 
the slope of the mean-camber line thus becomes 

■yfl^W [-Ai'-jz^eS (2A'+^/)] COS e+ 



(50) 



(51) 



VI— Af'S-^-'cosftfl 
The coefficients in this case are given by the rdations 



Thus for the airfoil in the tunnel the coefficient A^,' is a 
function of both the angle of attack and the shape of the 
camber line, but the functional relationship is altered from 
what it was in free air by the iadusion of terms proportional 
to <r. Furthermore, because of the appearance of the term 
involving Ao' in the second of> equations (51), the coefficient 
Ai is in this case also a function of the angle of attack, as 
well as of the camber-line shape. Since .^o' appears in this 
equation multiplied by the factor (7, the dependence upon the 
angle of attack is, however, secondary as compared with 
the dependence upon the shape of the camber hne. As in 
the case of the airfoil ia free air, the coefficients A„' for 
71^2 axe functions of the camber-line shape only. 

The chordwise lift distribution in the tunnel is given by 

dL' ^jdV 2dT' 



or in coefficient form. 



(52) 



In writing this equation the streamwise velocity gi-atlient 
which results from the wall-wake interference (equation 
(24)) is ignored. It can be shown that the inclusion of this 
variable would give rise to correction terms of the order trca'. 
Terms of this order are usually small as compared Avith the 
terms of order <r and tc/ considered in the theory and may 
therefore be neglected. 

It is apparent from equation (52) that, as in the case of 
the airfoil in free air, the lift distribution in the tunnel may 
be divided into two components. Now, however, the com- 
ponent which depends upon the angle of attack includes both 
the cotangent term and the first sine term. The component 
which is a fimction of the camber-line shape alone comprises 
the sine terms corresponding to 7i^2. Again, those two 
components could be denoted by the terms "additional" and 
"basic" in the sense previoxisly employed; however, smco 
the phrase "additional lift" already is so firmly established 
with reference to the distinctive cotangent term alone, this 
xisage does not appear advisable in the present cose. For 
this reason, the terms of the series will be referred to by 
reference to their form or their position in the series. 

Since it is the same airfoil which is being considered in 
both cases, equations (47) and (51) lead to the foUowng 
general relations between the coefficients in free air and in 
the tunnel: 

«_ vr=m4o= «'- vr=^A' + Yi=^(^' + ^^»') 



(63) 



A2=A3,' 

a^^a: 

In order to use these expressions to relate the characteris- 
tics of the airfoil m free air with those in the tunnel, it is 
necessary to choose some quantity or condition which Avill be 
maintained the some in both cases and relate the remaining 
quantities in accordance with this choice. If it were pos- 
sible, the ideal procedure would be to keep all the aerody- 
namic coefficients unaltered and to determine a correspond- 
ing relationship between the angle of attack in the tunnel and 
in free air. To do this it would be necessary to keep all 
pressure and frictional forces the same in both cases, which 
can be accomphshed only if the pressm-e distributions are 
identical. This would require that each of the coefficients 
An in equation (48) be equal to the corresponding coefficient 
A„' in equation (52). It is apparent from the second of 
equations (53), however, that this requirement cannot, in 
general, be satisfied. 

Although the pressure distribution cannot be maintained 
completely unaltered in the transfer from the tunnel to free 
air, the general relations (53) can be specialized in such a 
way that the essential character of the distribution is un- 
changed. It is apparent that the component of lift contribu- 
ted by the first, or cotangent, term ia equations (48) and 
(52) is different in form from that contributed by the scries 
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of sine terms. The cotangent component has an infinite 
value at the leading edge (6=0) and a relatively laorge chord- 
wise gradient of lift over most of the chord of the airfoil. The 
sine-series component is finite at all points and, for airfoils 
ordinarily encountered in practice, has a relatively small 
chordwise gradient, except possibly in the immediate vi- 
cinity of the leading or trailing edges. The cotangent 
component with its infinite peak pertains, of course, only 
to the hypothetical airfoil of infinitesimal thickness and 
zero leading-edge radius. For aU real airfoUs, the lift at the 
leading edge can never be infinite; however, even in this 
instance the lift distribution is characterized by a component 
the form of which is peaked near the leading edge and the 
magnitude of whicb varies markedly with the angle of attack. 
The magnitude of this component is a primaiy factor in 
determining the character of the pressxu-e distribution, and 
even a relatively small change in magnitude may cause con- 
siderable change in the TniniTmim pressure and in the chord- 
wise pressm-e gradients attained on the surface of the airfoil. 
Further, the aerodynamic chaiacteristics which depend upon 
these quantities, particularly the profile drag, TTinTiTrmTn lift, 
and critical compressibility speed, will be correspondingly 
altered. It follows that properly to correct airfoil data 
obtained in a wind tunnel to conditions in free air, the cor- 
rected quantities should correspond to the same magnitude 
of the peaked lift component as exists on the airfoil ia the 
tunnel. 

The requirement that the peaked component of Hft on 
the real airfoil shall be the same ia the tunnel and in free 
air can be expressed with reference to the assumed airfoil 
of infinitesimal thickness and camber by sotting equal 
to Ac' in equations (53). The first of these equations, which 
relates the angle of attack in the two conditions, then 
becomes 



(64) 



and the relations between the coefiidents are 

An=A;' 



Substitution from equations (56) into equation (48) 
gives 

P=4 [a' cot |e -jzJl? (2A'+^i') sin "+2 A' sin nfl] 



(55) 



or 



4(r 



P=P*- j^f^ (2A'+Ai') sm 



(56) 



Thus, if the angle of attack in the tunnel and the angle of 
attack in free air arc such as to satisfy equation (54), the 
84311ft— 50 12 



chordwise lift distributions will differ by an amount defined 
by the second term on the right-hand side of equation (56). 
The lift coe£Glcient for the airfoil in free air is 



which, after substitution from equations (44) and (48), can 
be integrated to give 



c,=ir (2^,+^,) 
The quai'ter-chord-moment coefficient is 



(57) 



=/>(H)^(f) 



which becomes after integration 



c« =— ;j(^i— 4s) 



(58) 



In usual wind-tunnel practice, the measured coefficients 
are referred to the apparent dynamic pressure g'. The lift 
distribution over the airfoil in the tunnel in terms of g' is 

p/=i dL'_2 £ 
2' dx Vqf dx 

Substitution from equation (49) gives 



P'=4 ^ j.4o' cot I e+YuAn' sin Tifl 



(69) 



The lift and moment coefficients of the airfoil in the tunnel 
as referred to the apparent dynamic pressure are then, 
respectively, 



c,'=J^'p'd(^)=^i {2A,'+A,') 



(60) 



and 



Relations between the coefficients in free air and in the 
tunnel can now be found with the aid of equations (55). 
Substitution of values from these equations into equation 
(57) gives 

c,=:r(2A'+4i') (l-IZ^) 

Substitution from equation (29) and neglect of correction 
terms of higher than the first order then give 

_ /Ji <r 2-(M0* . 

|i-iir^-[i_(M')«]s/* 

[2- (MOT [1-|-0.4(M0'] , 



166 



REPORT NO. 782 NATIONAL ADVISORY COMJVtlTTEB TOR ABRONAtmCS 



From equation (33) it can be seen that, to the first order, M 
may be replaced by M' in this equation. The final equation 
for the correction of the measured lift coefficient is therefore 

_ ,U__o_ 2-{M'y . _ 



[2- (MOT [1+0.4(M0^ 
1-(M')* 



rc/J 



(62) 



Similarly, substitution of values from equations (55) into 
equation (58) gives 

1 U.'-Ai')+\ (2A'+^i') 

To the order of approximation previously employed, the final 
equation for the correction of the measured moment coeffi- 
cient can be written 



2-(M'y 



[2-(M0»][l+0.4(M01 

i-(,M'y 



A<r- 



'^'^^ 4[1-(M')T 



(63) 



The corresponding angle of attack in free air can be found 
from equation (54). Combination of equations (60) and 
(61) gives 

^ ^{c,'+4c. 



2a- 



To the first order, equation (54) then gives for the corrected 
angle of attack in radian measiu'e 



2W1-(M')' 



4 



or in degrees 



, , 57.3(r f , 



(64a) 
(64b) 



Numerical values of the compressibility factors appearing in 
equations (62), (63), and (64) are given in table II. 

It shoiild not be implied from equations (64) that the 
general iaclination of the stream at the position of the airfoil 
in the tunnel is actually different from what it would be if the 
walls were not present. The equations indicate rather that, 
with regard to the magnitude of the cotangent component of 
lift distribution, an airfoil at a given angle of attack in the 
tunnel behaves as though it were at a different angle in free 
air. This difference occurs because the tunnel walk give rise 
effectively to a ciiange in the curvature of the stream at the 
position of the airfoil. 

As was indicated previously, the essential character of the 
pressure distribution over a given airfoil wiU be the same 
in tihe tunnel and in free air, provided the magnitude of 



the cotangent lift component is the same in both coses; that 
is, provided the angles of attack are such as to satisfy 
equations (64). The exact shape of the pressure distribu- 
tions, however, will still differ slightiy for two reasons: 
(a) The interference between the lift and the tunnel walls 
causes a difference in chordwise lift distribution as required 
by equation (56), and (b) the interference between the wake 
and the walls gives rise to a longitudinal pressure gradient 
defined by equation (25). The effect of tiiese two influences 
upon the remaining airfoil characteristic, the profile-drag 
coefficient, must be considered. 

As given by equation (56), the chordwise lift distributions 
in the tunnel and in free air differ by an amount 

AP=P*-P=j3^ (2^'-|-4i') sin B 

which, by virtue of equation (60), may be written to the 
first order as 



» n 4 <r I • . 



(65) 



The changes in peak pressure and pressure gradient brought 
about by this increment of lift distribution, unlike the 
changes which would accompany even a minor alteration of 
the cotangent lift component, are ordinarily small. At low 
Mach nxunbers the drag depends primarily upon the char- 
acter of the flow in the boundary layer, and, since this flow 
will not ordinarily be altered greatiy by these small changes 
in the pressure distribution, the increment of lift distribution 
should have only a small effect upon the profile drag. At 
high Mach numbers the dtag is determined primarily by the 
total-head losses in the shock waves which appear after the 
critical Mach number is reached; that is, after the local 
speed of sound is obtained at the minimum pressm-e point on 
tie airfoil. The critical Mach nxmiber is usually reduced by 
the change in peak pressure accompanying the change AP 
in lift distribution, but it can be shown that this reduction 
is ordinarily very small. It is reasonable to expect that the 
change in profile drag at a given supercritical Mach number 
is correspondingly small. These changes are discussed in 
fiirther detail later in the report, but for the present it may 
be assiuaed that the difference in chord^vise lift distribution 
between the tunnel and free air has only a negligible effect 
upon the profile drag. 

For usual airfoils and drag coefficients, the longitudinal 
pressure gradient defined by equation (25) is also small, and 
its effect upon the boimdary-layer flow and hence upon the 
frictional drag of the airfoil may be neglected. It wiU, 
however, increase the pressure drag by an amount which is 
comparable to differences already retained in the corrections 
to the lift and moment. This increase va. pressiu-e drag 
must be subtracted from the drag measm-ed in the tunnel to 
obtain the true profile drag of the airfoil in free air. 

Glauert has shown (reference 2, pp. 62-63) that in an 
incompressible fluid the drag experienced by an airfoil as 
the result of a streamwise pressure gradient is, in the nota- 
tion of this paper. 



tJj= s A<r -f-— Ao- -f- 

8 ax T ax 



(66) 
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In appendix A, it is shown that to the order of the linear 
theory the relationship between drag and streamwise pres- 
sure gradient is unchanged by the effect of compressibility. 
However, the use of factor A in equation (66) implies, as in 
the case of the velocity correction due to thickness, a second 
order correction which should properly be modified for the 
effect of compressibility. This modification which has been 
determined by Ludwieg, reference (10), is not included here 
for the reason previoxisly cited. Substitution in equation 
(66) of djpldx from equation (25) then gives for the drag due 
to the interference between the wake and the walls 



ATI I l+OAQd'Y 

coj—Ci g;c \ pT^]^?]^ 



A<r 



The true profile drag of the airfoil in free air is then 
D=D'-6D 



=Ci'q_'c 



, 14-0.4(MT , 
^ [1-(M')T'' 



and the corresponding drag coefficient referred to the true 
dynamic pressure is 

Substitution from equation (29) gives for the final correction 
to the measured drag coefficient 

//, 2-(M0' . 1+0.4(M0' , 



[2-(M0^ [\+QA{M')^ 
1-(M')* 



(67) 



It will be noted that, of the two correction terms involving 
A<r in this equation, the first appears as a resiilt of the change 
in dynamic pressure occasioned by the interference between 
the walls and the airfoil thickness; the second represents the 
effect of the pressure gradient induced by the interference 
between the waUs and the wake. The correction term con- 
taining TCi appears as a result of the change in dynamic 
pressure caused by the wall-wake interference. Numerical 
values of the functions of M' which appear in equation (67) 
are given in table n. The corrected drag coefficient cor- 
responds, of course, to the corrected lift and moment co- 
afficients as given by equations (62) and (63) and to the 
corrected angle of attack as given by equation (64a) or (64b). 

The drag correction of equation (67) was detennined 
particularly for drags measured with a balance and, as de- 
rived, is not necessarily correct for drags measured by the 
wake-survey method. It can be shown, however, from 
theoretical considerations of momentum and continuity in 
a two-dimensional-flow tunnel that for normal ratios of 
airfoil chord to tunnel height, the ordinary ts^pe of wake 
smrey derived for free-air conditions gives, when apphed 
in the tunnel, a value of the drag equal to that measured 
by the balance except for a neghgible difference of less than 
one-half of one percent. Equation (67) may thus also be 



used to correct the drag coefficients determined by the wake- 
survey method. 

It should be noted that no correction to the dr^ has been 
made for any pressure gradient which may exist inherently 
in the tunnel as a result of the streamwise growth of the 
boundary layer on the tunnel walls. Most modern tunnels 
are constructed so that this pressure gradient is sensibly 
zero; however, if such a gradient does exist and its magnitude 
is known, an approximate correction to the airfoil drag can 
be made by means of equation (66). 

There remains the necessity for correcting the measured 
pressure distribution over the siuface of the airfoil. The 
pressiure at any point on the airfoil is conveniently expressed 
by the pressure coefficient 8i defined by 



or by the pressure coefficient P| defined by 



(68) 



(69) 



where 'pt is the local static pressure on the surface of the air- 
foil and H, and 2 ox&, respectively, the total head, static 
pressm-e, and dynamic pressure of the undisturbed stream. 
As indicated in reference 14, in a compressible stream, 



ff=3)+2(l+,) 
where (l-H'?) for air (7=1.4) is defined by the series 

,M* Af 



l-|-r?=l- 



4 ^40 ^1600 



4- 



(70) 



(71) 



M being the Mach number of the stream. lYom these rela- 
tions it is readily shown that 



(72) 



A curve of (I-I-17) versus M, as calculated from equation 
(71), is given in figure 4. 

In reference 6 a method is presented for the determination 
of the pressure distribution around an airfoil in an incom- 
pressible stream when the lift distribution along the chord 
and the pressure distribution over the base profile are known. 
The upper- and lower-surface pressures at any chordwise 
station x are given ia. coefficient form by 



P ~2s. 




[(i-p.)+ipj 
_ [(l-p.)-^pj 



(73) 



where Py is the pressure coefficient on the base profile at x, 
and P is the coefficient of lift per unit of chord at x. By 
following the basic reasoning of reference 6 and assuming 
that the induced velocities at the siurface of the airfoil are 
small as compared with the velocity of the undistm-bed 
stream, it is readily shown that equations (73) may also be 
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applied to the pressure distribution in a compressible stream. 

In such, application, the values of Pp, P^, Pf, and P must all 
correspond, of course, to the same free-stream Mach number. 

The measured pressure distribution is now readily cor- 
rected for tbe effect of the tmmel walls. It is only necessary 
to refer the measured pressure coefficients to the true instead 
of the apparent dynamic pressure and remove the effect of 
the Uft distribution represented by equation (65). Strictly 
speaking, correction should also be made for the pressure 
gradient due to the waU-wake interference; however, in 
practical tests sucb correction is smaU and may be neglected. 
The detailed procedure is then as foUows: 

(1) The apparent upper- and lower-smrface pressure 
coefficients 



and Srf-- 



are obtained from the experunental results lor the various 
chordwise stations. 

(2) These pressm:e coefficients are referred to the true 
dynamic pressure by means of the equations 



\-{M'y 



12- (M')^ [1+0-4(MO'] 



(74) 



(3) Tbe quantities (1 — Pu*) and (1 — P^*) are determined 
in accordance with equation (72) as 



l-Pu*=Sv*-v] 
l-Pi.*=Sr.*-v\ 



(75) 



where 17 is determined by figiu^e 4 for the true Mach munbor 
as given by equation (33). 

(4) The chordwise lift distribution iu the tunnel is found 
from 



(76) 



(5) The chordwise lift distribution iu free air is deter 
mined from equation (65), which may be written 



P=P* 



'l-(M')* ' ' 
where P, is given by 



(77) 



(78) 



This quantity, which is termed the "interference lift distri- 
bution," is seen to be elliptic in form. Values of P, at stand- 
ard chordwise stations are given in table III. 

(6) The quantity (l—P/), where P/ is the base-profile 
pressure coefficient in free air, is given by the equation 



(70) 



which is obtained by combining equations (73) . 

(7) The values of P and (l—P/) being known, the uppor- 
and lower-surface pressure coefficients Pv and Pz, are deter- 
mined from equations (73). If desired, the corresponding 
coefficients Sa and Sl can be found from equation (72) . 

The corrected pressure distribution obtained by this 
method corresponds to the corrected angle of attack as given 
by equation (64a) or (64b) and to the corrected lift and 
moment coefficients as given by equations (62) and (63). 

It has been mentioned previously that the correction to 
the angle of attack appearing in equations (64) does not 
represent an actual rotation of the stream direction. This 
fact is rmphcit in the derivation of the equations, but it can 
also be demonstrated by simple considerations of force and 
momentum. For this purpose it is sufficient to consider a 
simple incompressible potential flow in the tunnel and ignore 
the effect of the profile drag. Assmne for the time being 
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that, because of the interference between the airfoil and the 
tunnel walls, the general direction of the stream at the air- 
foE is inclined from its original direction parallel to the 
tunnel walls. For potential flow the resultant force acting 
on the airfoil must be at right angles to the local direction 
of the stream. The airfoil thus would be acted upon under 
the assxuned conditions by a component of force parallel to 
the center liae of the tunnel and would in reaction exert an 
equal and opposite force on the flow. Since the timnel walls 
cannot in a potential flow exert a force parallel to the center 
line, this longitudinal force would have to be balanced by a 
difference of pressure or momentum between two stations 
in the tunnel, one upstream and one downstream from the 
airfoil. If the stations ore taken far enough from the air- 
foil that its induced velocities are negli^ble, conditions 
across the timnel are uniform at each station. It then fol- 
lows from considerations of continuity of the incompressible 
flow in the tumel that the conditions at the two stations are 
identical, and no difference of pressure or momentum is 
possible. Thus the original assvunption of a general rota- 
tion of the stream direction at the position of the airfoil is 
untenable. This conclusion is not changed by the effects of 
fluid compressibility. Furthermore, the fact that the in- 
troduction of the profile drag and the accompanying wake 
causes a pressure difference between the two stations like 
wise does not alter the result, as the wake effects are con- 
sidered in the theory to be superposed on the potential-flow 
field. Thus, the angle correction appearing in equations 
(64) must be due to some cause other than a general inclina- 
tion of the stream. As previously pointed out, it is actually 
due to an effective change in the curvature of the sttream at 
the position of the airfoil and is a direct consequence of the 
requirement that the airfoil in this stream shall have the 
same cotangent component of lift distribution as does the 
airfoil in free air. These considerations are important in 
the proper interpretation of drag measurements from a two- 
dimensional-flow tunnel. 

In the development of the correction to the measured drag 
coefficient, it was assumed that the increment AP in chord- 
wise lift distribution between the tunnel and free air has 
only a negligible effect upon the profile drag. A better idea 
of the nature of the effect can be had by further examination 
of the difference between the two cases. It follows from 
equations (65) that, if the angles of attack in the tinmel and 
in free air are related as required by equation (54) or (64), 
the transposition of a given airfoil from free air to the tmmel 
is equivalent to increasing the coefficient Ai for the airfoil 
in free air by an amount 

which can be written to the first order as 

As can be seen from equation (46), this can be accomplished 
by maintaining the angle of attack imaltered in free air and 



changing the ordinate of the mean-camber hne at every 
point by an amount Ay^ such that 



(£(Ay,) ^l a_ 



Ci cos B 



(80) 



The value of Ay^ as a fraction of the chord is then 



t-f-^<f)+« 



which after substitution from equations (44) and (80) can 
be integrated to obtain 



Ay,^ 1_ 0 

c Sir Vl— M* 



c, cos 2d+G 



The constant of integration G is determined by the condition 
that At/Jc=0 at 9=0 and 6=%. The equation for the 
change in the camber line then becomes finally 

^•=^7T^4(f)-(D"] 

This is the equation of a parabola with vertex at the mid- 
chord point and has the same form as the equation for the 
camber line of an NACA conventional airfoil with maximum 
camber at the midchord point (reference 16). The maximum 
change in camber is 



Ayc\ _1_ 

€/ /max 



(82) 



Th\is, if the angles of attack of the airfoil in the tunnel and 
free air are adjusted as required by equations (64), the wall 
interference in the tunnel has the same effect upon the chord- 
wise lift distribution as would an increase in camber in free 
air. 

As a possible instance of a test for the determination of the 
drag of an airfoil of large chord at a low Mach number and 
low lift coefficient, consider the case of an airfoil in a tunnel 
providing a chord-height ratio of 0.5. The value of o- is 
then 0.051. Assume that the angle of attack a' in the tunnel 
is adjusted as required by equations (64) to correspond to 
an angle a giving a lift coefficient c, of 0.3 in free air. As- 
suming that the Mach niunber is siifficientiy low that the 
effect of compressibiUty may be neglected in computing the 
tunnel-wall jcorrections, the change of maximum camber re- 
quired in free air to duphcate the effect of the tmmel walls is 
given by equation (82) as 



( 



^ /nuuE 



^ (0.051) (0.3) =0.0012 



An estimate based upon experimental data has been made of 
the effect upon the profile drag of a change in camber of 
this magnitude for an NACA conventional airfoil of moderate 
camber and 16-percent thickness with maximum camber at 
the midchord point. The result indicates that neglecting 
the effect upon the proffie drag of the change in lift distribu- 
tion caused by the tunnel walls leads in this case to an error 
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in the final corrected drag coefficient of less than 0.0001. 
This is within the usual limits of experimental accuraxrjr. 
The correction terms included in equation (67) amount in 
this instance to approximately 0.0004. If the chord-height 
ratio were increased to 1.0, the error in the drag coefficient 
would he increased to 0.0004, which is well outside the limits 
of experimental accuracy. This indicates the desirability 
of limiting the chord-height ratio if accurate measurements 
of the profile drag are desired, even at low values of the lift 
coefficient and Mach nxunber. At higher values of the lift 
coefficient or Mach number the permissible chord-height 
ratio must be reduced correspondingly. 

The foregoing comparison is based upon the specific case 
of an airfoil with maxumun camber originally at the mid- 
chord point and is not necessarily applicable to other types 
of airfoils. For families of airfoils which have a smaller 
variation of drag with camber than do the NACA conven- 
tional sections, tiie error introduced by neglecting the effect 
of the change in hft distribution is correspondingly less. 
In any event, satisfactory accuracy can be obtained in the 
measurement of drag at low lift coefficients and Mach num- 
bers by keeping the chord-height ratio within a suitable 
limit — say 0.7. A possible exception is an airfoil having an 
essentially flat pressure distribution in the region of tran- 
sition from laminar to turbulent flow in the boundary layer. 
In such a case the changes in pressure gradient may shift 
the point of transition and noticeably alter the profile drag; 
however, for any such sensitive airfoil, alterations from this 
source are no more serious than similar changes which may 
accompany the small variations in pressiure distribution 
caused in any practical application by irregularities in 
construction. 

Some measure of the effect of the increment AP ia chord- 
wise lift distribution upon measurements of airfoil charac- 
teristics at high Mach numbers can be obtained by calcu- 
lating the change in critical Mach number caused by this 
increment. Such a calculation has been made for an airfoil 
with miniTmim pressure originally at the midchord point. 
Since the increment AP is a maximum at midchord, this 
represents the worst possible case as regards the change in 
critical Mach ninnber. For a chord-height ratio of 0.25, 
which is considerably laiger than that ordinarily employed 
in tests at high Mach numbers, the critical Mach number 
was foimd to be reduced by approximately 0.001 at a lift 
coefficient of 0.3. A change of this magnitude is insignifi- 
cant. It may be expected that the accompanying change 
in the aerodynamic coefficients in the vidnilry of the critical 
Mach number wiH be correspondingly small. 

THE PHENOMENON OF CHOKING 

CJonsider the compressible adiabatic flow of a fluid in an 

elementary stream tube of varying area Ai, as shown in 
figure 5(a). Continuity of flow requires that the product 
PiViAt be constant, where pi, Vi, and At are the local values 
of density, velocity, and area, respectively, at any station 
along the tube. In consequence, the logarithmic derivative 
must vanish; that h, 



dpi,dV(.dAi - 



(83) 



A, 



Pt 



(a) 




FiOTTBE S.— Velodty dlstribntlon In an elementary stream tubo. 

Bernoulli's equation for compressible flow requires that 



dpt^ 



Pi 



V,dVt 



(84) 



where pi ia the local pressure. Defining Vc, as the local 
velocity of soimd, then 



dpi 



V » 



so equation (84) becomes, after substituting the value of 
dpi in that equation, 

dp,_ V,dV, ^,dV, 

77 V:? -^'TT 

where M| is the local Mach nxmiber. 
Substituting this relation into equation (83) gives 



il~M,^ -y^- ^ 



(86) 



From this weU-known relation it is seen that at subsonic 
speeds the usual behavior associated incompressible 
flow is obtained; namely, that as the area increases the ve- 
locity decreases. At supersonic speeds, however, the be- 
havior is reversed in that as the area increases the velocity 
increases. When the local Mach number is unity it is seen 
that dA=0; that is, if the velocity of sound is attained in the 
tube it can only be attained where the area has its minimum 
value. 

When the local velocity of sound is attained at the mini- 

vamn. area section, the local Mach nmnber at any other sec- 
tion, determined by the ratio of the area at that section to 
the minimum area, may be less or, in some cases, greater 
than unity depending upon the conditions promoting the 
flow in the tube. The nature of such flows can be studied by 
considering the change in the character of flow in the stream 
tube of figiu-e 5(a) as the downstream pressm-e is do- 
creased with respect to the upstream pressiu-e pi. If 
Pi — 3?2 is small so that completely subsonic flow is maintained 
in the tube, the nature of the velocity variation along the 
tube is that usually associated with incompressible flow as 
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seen in curve I of figure 5(b). When pi — p^ is increased so 
that sonic speed is just reached in the TninimuTn area section , 
the variation of velocity along the tube becomes that shown 
in curve II. Any further decrease of the pressure cannot 
alter the flow upstream of the minimum area, since the velo- 
city at the minimum section cannot exceed the velocity of 
soimd. The only effect of decreasing the downstream pres- 
sure is to promote a supersonic flow region downstream of the 
niinimiim area, as shown by curve III of figure 5(b). This 
region is terminated by an abrupt return, through a com- 
pression shock wave, to subsonic flow. The position of this 
terminal shock wave must be such as to bring about the 
necessary conversion of kinetic to thermal enei^ that is 
required to promote the downstream pressure p^. For pres- 
ent purposes, the most important point concerning the flow 
as described is that when the velocity of sound is attained 
at the Tnim'Tmim area section, no fm-ther increase in the 
flow rate can be made regardless of the extent of the super- 
sonic flow region downstream of this section. When this 
maximimi flow rate has been reached, the stream tube is 
said to be "choked." 

What has been said concerning the choking of a single 
stream tube applies to the complete system of stream tubes 
past an airfoil mounted in a two-dimensional-flow tunnel, 
as shown in figure 6. That is to say, when the velocity of 



M'c), 




FiouBE 6.— Lines of sonlo speed at ths position of tbe airfoil after choking 

the undisturbed flow far upstream in the tunnel reaches a 
certain value, sonic velocity is attained at the point of mini- 
mima area of each elementary stream tube betwe^ the airfoil 
and the upper wall of the tunnel. It is important to note 
that the locus of the points of minimum area of the separate 
stream tubes does not necessarily coincide with the shortest 
line between the airfoil and the upper wall. This is illus- 
trated in figure 6, where the line A represents the shortest 
distance between the airfoil and the wall. If the conditions 
of flow were uniform across the stream at each chordwise 
station, the flow between the airfoil and the wall would be 
the same as in a single elementary stream tube, and sonic 
velocity would necessarfly be attained along line A. In the 
actual case, however, the flow is two-dimensional, and sonic 
velocity is attained along some hne, such as line B, not coin- 
cident with A. A sitrular situation exists in the space 
between the airfoil and the lower wall of the tunnel, where 
the sonic velocity is attained along some hne D. As before, 
this lino does not necessarily coincide with line 0, the short- 
est hne which can be drawn from the lower surface of the 
airfoil to the lower wall. (In order to avoid an apparent 
contradiction with the requirements of continuity, it must 
be kept in mind that the velocity vector is not, in general, 
perpendicular to either hues A and C or B and D.) Sonic 



speed is generally not attained coincidentaUy along lines B 
and T>. Once it is attained along both these lines, however, 
the rate of flow past the airfoil ia the tunnel can undergo no 
fiurther increase. The Mach number of the flow ahead of 
the airfoil then has its maximum attainable value. This 
value is described as the "apparent choking Mach mmiber." 

In practice, the lines of sonic speed lie very dose to the 
lines defining the shortest distance between the airfoil and 
the tunnel walls. For purposes of analysis, it will be as- 
sumed that they are coincident, that is, that Unes B and D 
coincide," respectively, with lines A and C. Under these 
conditions, the calculated rate of flow in the tunnel (which 
must in any event be equal to the rate of flow across lines 
A and C) vnR be somewhat greater than that which actually 
exists when the lines B and D have their true positions. 
The assumption of imidimensional flow will thus lead to a 
computed choking Mach number, which is slightly greater 
than the theoretically correct value. 

On the basis of the foregoing assumption a relationship 
between the model size and the chokmg Mach mmiber can 
be obtained from elementary considerations. The velocity 
V and density p' of the flow far forward of the model, where 
the cross-sectional area is A', are constant across the stream. 
The velocity V„ and density p„ across the sonic-speed lines 
B and D of figure 6, where the area has the minimum value 
A„, are again constant across the stream. The velocity 
V„ is the local sonic speed Ve^ so that the equation of 
continuity becomes 

Assuming adiabatic relations, the density and velocity 
terms can be related to the Mach number far upstream, 
which is now the apparent choking Mach number. The end 
result is that the ratio of the area of the undisturbed stream 
to the TniniTrmTTi flow area can be expressed -in terms of the 
apparent choking Mach number M'eb as 



1 



1 + 



The area ratio is clearly 
A 



7+1 



-1] sTFi) 



(86) 



A„ h—tj, 



where h is the tmmel height and the projected thickness 
of the airfoil normal to the flow direction. For reasons 
which will be evident later, the projected thickness in t.Viia 
relation will be replaced by an "effective" thickness t,. 

Taking the value for 7 for air as 1.4, equation (86) 
becomes 



(87) 



In figure 7, the ratio tjh is plotted as a function of the 
apparent chokuig Mach number. The region above the 
curve represents an impossible state of flow. As a matter 
of interest the residts are shown for the supersonic- as well 
as the subsonic-flow regitne, although for the purpose of this 
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FlouBE 7.— Choking Mach number as a fanatlon of effective-thlskness to tannd-helght ratio. 



report only the subsonic choking Mach niunbers "will be 
considered. 

In writing equation (87), the projected thickness was 
replaced by an effective thickness. If choking occurred as 
was assumed in the preceding analysis, then the effective 
thickness determimng choking would be, of course, the pro- 
jected thickness. In any real case, although the effective 
thickness may never be less than the projected thickness, it 
may be greater for two reasons. First, if the angle of attack 
is sufficiently large in absolute value, one of the lines B or 
D may move downstream of the trailing edge because of the 
continued contraction aft of the airfoil of the portion of the 
stream passing that line. Second, since on any aerody- 
namic body there exists, because of the action of viscosity, 
a boimdary layer wherein the velocity must be reduced 
below the velocity in. the otherwise imaffected flow field, 
it follows that the velocity of soimd cannot be attained at 
those points close to the airfoil surface on the lines B and 
D of figure 6. 

To estimate the choking Mach number in any practical 
case, it is necessary to assume that the effective thickness 
is equal to the projected thickness of the airfoil. Because 
of the possible contraction of part of the streajgi aft of the 
airfoil, as weU as of the assumption that unidimensional 
flow exists as previoxisly described, this procedure will lead 



to a computed choking Mach number which is greater than 
the theoretically correct values for an ideal, incompressible 
fluid. Further, the influence of the boimdary layer will 
cause the actual choking Mach number to be even less tlian 
this theoretically correct value. Thus the use of the pro- 
jected thickness in the computation may be expected to lead 
to an overestimation of the choking Mach number. 

The effect of the boimdary layer in this regard may best 
be illustrated by the case of a flat plate set at zero angle of 
attack in a two-dimensional-flow wind tvmnel. Since the 
projected thickness is zero, the previously developed theory 
would indicate that no choking would occurr in this case. 
Actually, because of the fact that the plate has a boimdary 
layer and an accompanying wake, choking does occur, as is 
shown in the following discussion. 

It was seen in the section on wake effect, wherein the 
effect of confining the wake of a body experiencing drag was 
considered, that when the influence of the wake spreads to 
the waUs so that a imiform velocity field again exists, the 
temperature at this downstream position is related to the 
temperature upstream of the model by equation (13). Using 
equation (9), tiie ratio of the corresponding velocities may bo 
seen to be 



V'- (7+l)(M')" 



(T+1)(MT 



(88) 



The velocity ratio is imaginary when the sign of the group 
of terms imder the radical is negative. The fimctional 
relationship beween the choked Mach number and the drag- 
density factor TCa, found by equating the terms under the 
radical to zero and solving the resulting equation, is thus 
determined as 



TCa' 



(89) 



where, as before, t=^ (j^- Setting 7=1.4 for air gives 

^ ,_1±1A(M^ I r l-(M-^ )^~Tl , . 
^''^ - 2.S(M'^y P-V^ Ll+1.4(M'^)«J 1 

A graph of this function is shown on figure 8. The effect of 
drag on choking for supersonic as well as subsonic wind 
timnels is shown as a matter of interest. 

The manner in which drag promotes choking may be com- 
prehended by examining the variations of the ratios V" jV 
and T"IT' in equations (88) and (13) as the value of tc/ is 
increased. In the case of the subsonic ■wind tunnel, the 
effect of increasing rCa is to increase V'lV. On the other 
hand, T"IT' and hence Vc'IVc are reduced. Consequently 
M"IM' is increased. In the case of the supersonic wind 
tunnel the effect of increasing tc/ is to decrease V'lV and to 
mcrease T"IT' and hence Vc'^Vo- Consequentiy, M" jW 
is reduced ui this case. In both cases choking occurs when 
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FiouBE 8.— OhoUng Iklaoh nomber as a fonotian of drag-density botor. 

the value of TCi is such as to make fihe downstream Macih 
number M" unity. 
There is one definite limitation of the previous analysis in 

that it was assumed that the effective tunnel area remained 
constant at least imtil the wake had spread to the walls so 
that uniform flow conditions were obtained across the stream. 
Such a condition does not prevail in any conventional wind 
tunnel, nevertheless the results are useful in providing ap- 
proximate values for the effect of drag as it determines 
choking. For example, a flat plate having an apparent drag 
coefBcient of 0.007, if the chord-height ratio were 0.5, 
would choke a subsonic wind tuimel at a Mach number of 
0,95 if choking occurred as assmned in the analysis. The 
serious influence of drag on choking for airfoils for which the 
drag coefBcient may be-many times this value is evident. 

To summarize, it has been shown that choking will occur 
in a wind tunnel as a result of the confinement of the flow 
caused by the presence of the model and its wake. In the 
case of airfoils of normal thickness, choking will usually 
be determined by the effective dimensions of the body — that 
is, by the actual dimensions modified for the effects of 
boundary layer and stream contraction aft of the airfoil as 
previously described. Properly, the boimdary-layer effect 
is a drag influence, but since its contribution is usually small 
it is most convenient to classify such confinement effects 
along with those due to the physical airfoil dimensions. In 



the case of very thin airfoils, at small angles of attack, 
choking wiU usually result from the confining effect of the 
wake rather than the effect of the airfoil thickness. 

Once the choking Mach number is reached, no further in- 
crease in tunnel power can affect the apparent Mach number. 
Such an increase will only serve to extend the supersonic flow 
region downstream of the lines of sonic speed. The forces 
experienced by the airfoil at choking thus vary depending 
on the power input to the wind tunnel. 

As a final consideration it should be noted that the flow in 
the timnel at choking does not correspond to any real flow 
over an airfoil in free air. Since the choking Mach niunber 
approaches unity as the tunnel height h becomes infinite, 
flow in the tunnel at choking, if it is to correspond to any flow 
in free air, must correspond to the flow that would occm- 
around an airfofl in a free stream moving at the velocity 
of sound. It can be demonstrated, however, that such 
a correspondence is impossible. Experimental evidence 
indicates that the flow conditions existing in the tunnel at 
choking are essentially steady state. That the flow about 
an airfoil in a free stream having the velocity of sound 
cannot be a steady-state flow can be readily shown. For 
instance, it was demonstrated previously that in any stream 
tube the velocity of soimd, if it is attained at all, must be 
attained at the minimum area section. That is to say, the 
rate of flow per unit area is a maximum where the velocity is 
the velocity of sound. Now, presuppose a steady-state 
flow in the stream tubes in the vicinity of an airfoil in free air 
when the stream velocity is sonic speed. If the velocity 
either increases or decreases as the flow passes the airfoil, the 
stream tubes must expand. This is clearly impossible, since 
the disturbance to the flow would then increase continuously 
as the distance from the airfoil increases. "On the other 
hand, if the velocity remains the velocity of sound in each 
stream tube, the streamlines wiU then have the same shape 
at all distances from the airfoil. Also, the pressure wiU 
remain constant throughout the entire flow fidd. This is, 
of course, impossible, since pressure differences are necessary 
to promote the required changes in the direction of flow 
past the airfoil. A steady-state flow similar to that observed 
in the tunnel at choking therefore cannot exist in free air 
at a free-stream Mach munber of imity. Thus at the choking 
Mach number, the flow at the airfofl. in the tunnel cannot 
correspond to any flow in free air. It f oUows that, at choking, 
the influence of the tunnel walls cannot be corrected for. 
Further, in the range of Mach numbers close to choking, 
where the flow is influenced to any extent by the incipient 
choking restriction, any correction for wall interference must 
be of doubtful validity. 

That the flow at or dose to choking cannot be corrected 
for the interference effects of the tunnel walls may be rea- 
soned from another point of view. The assmnption that it is 
permissible to correct wind-tunnel test data for the influence 
of the walls is justified only when the influence on the flow 
near the model is of such a uniform nature as not to alter 
the general character of the flow materially from some 
corresponding flow in free air. For instance, a velocity 
correction for wall interference may be apphed with con- 
fidence only if the velocity increment resulting from such 
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interference is constant or nearly constant over that portion 
of the flow field wherein the influence of the model on the 
flow is important. Viewed in this light, it is clear that at or 
close to choking no correction can properly be applied, since 
an important influence of the model on the flow is felt over 
a range extending dose or up to the walls, within which 
range the influence of the walls on the flow is not at all 
uniform. 

It is thus clear that the equations which have been derived 
for correcting the test data obtained in a subsonic two- 
dimensional-flow wind tunnel for the effects of wall inter- 
ference, cannot apply at the choking Mach number nor for 
a range of Mach numbers below the choking value. More- 
over, when the model is not symmetrically disposed, the flow 
wiU, in general, attain sonic velocity across the stream on one 
side of the airfoil before it does on the other. In such cases, 
it is to be expected that the range of Mach nvmibers below 
choking for which the corrections are invalid is extended 
over that which would occur with a more nearly symmetrical 
flow pattern. 

DISCUSSION 

There is, at present, only a very limited amoimt of experi- 
mental data available which can be considered satisfactory 

for determining the accuracy of the theoretical interference 
corrections derived in this report. Moreover, none of the 
available data were obtained at sufficiently high Mach 
numbers to permit an evaluation of the accuracy of the 
theory with regard to the effect of compressibility. 

In figmre 9 are shown the experimentally determined 
variations of lift coefficient with angle of attack for several 
NAOA 0012 airfoils having different chord-height ratios. 
The data for those models for which the chord-height ratios 
are 0.25, 0.5, and 0.8 were obtained from tests in the 7-by 
10-foot wind tunnel at the Ames Aeronautical Laboratory. 
These models were of 6-foot span moimted across the 7-foot 
dimension of the test section; 6-iach-span dimimy ends were 
used in an attempt to obtain two-dimensional flow. A gap 
of about K» incli occurred between the test panel which was 
connected to the balance frame, and the dmnmy ends, which 
were fastened to the timnel walls. The lift was detenmned 
both from force tests and by integration of chordwise pres- 
sure distributions at a section close to midspan. The data 
presented here are those obtained from the pressure distri- 
butions. The data for the model for which the chord-height 
ratio is 1.0 were obtained from tests in the low-turbulence 
wind tunnel of the Langley Memorial Aeronautical Labora- 
tory. This 3-foot-span model was fastened directiy to the 
side walls of the tunnel such that no air gap existed, and the 
ift was determined from measmrements of the reaction on 
the roof and floor of the test section. The test results for 
the various models are shown, uncorrected for tunnel-wall 
interference, in figiu-e 9(a). In figure 9(b), are shown the 
same data corrected for wall interference by means of equa- 
tions (62) and (64b). For aU the models, the correction 
term depending upon tCj' is negligibly small. The test 
Reynolds numbers range from 2,000,000 to 6,000,000. It is 
seen that the corrected data obtained with the models for 
which (c/h) equals 0.25, 0.5, and 1.0 agree well with one 
another and with the section lift characteristics as obtained 
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from tests in the NACA variable-density wind timnel (refer- 
ence 15). The data obtained with the model for which (c/A) 
equals 0.8, when corrected, indicate a lower lift-cuxve slope 
than do the other data. This is thought to be due to the 
effect of air leakage through the gaps at the ends of the test 
panel, the influence of which may be expected to become 
more pronounced as the chord of the airfoil is increased 
relative to the span. 

In this regard, unreported tests in the Langley low- 
turbulence wind tunnel have shown that the presence of any 
gap through which leakage can occur will influence the aero- 
djTiamic characteristics to a surprisingly marked extent. 
This fact was also demonstrated by the Ames Laboratory 
tests on the NACA 0012 airfoils. A comparison of the lift 
characteristics obtained from balance measurements with 
those derived by integration of the pressure distributions, 
which are those given in figure 9, showed the lift-curve slopes 
for the former to be definitely lower than those for the latter. 
This indicates that the lift near the center of the test panel 
exceeded that at the sections near the gaps; that is, that the 
flow was definitely not two-dimensional. 

In figiu-e 10(a) is shown the experimental variation of 
lift coeflicient with angle of attack for an NACA 23012 for 
which (c/A) equals 1.0. These data were obtained in the 
Langley low-turbulence wind tunnel at test Reynolds num- 
bers of 4,560,000 and 6,450,000. The same data corrected 
for tunnel-wall interference by means of equations (62) and 
(64b) are shown in figure 10(b), together mth section lift 
characteristics as obtained in the variable-density -wind tun- 
nel at an effective Reynolds number of 5,000,000 (reference 
16). The corrected data are seen to be in excellent agree- 
ment with the results from the variable-density timnel. 

In figure 11(a) is shown the variation of quarter-chord- 
moment coefficient with lift coefficient for the NACA 0012 
airfoils as obtained from the 7- by 10-foot wind-tunnel tests 
previously described. In figure 11(b) are shown the some 
data as corrected for the interference of the tunnel walls by 
means of equations (62) and (63). The section moment 
characteristics for this airfoil as obtained from tests in the 
variable-density wind tunnel (reference 15) are also shown 
for comparison. It is seen that the corrected data are in 
fair agreement with the data from the variable-density wind 
tunnel, except for the model for which (c/h) equals 0.8. 
It is believed that this disagreement is again due to the effects 
of air leakage throiigh the gaps between the test panel and 
the dummy ends, and not to any shortcoming in the theory. 

In figures 12 and 13, the uncorrected and corrected projBIe- 
drag coefficients for six symmetrical bodies at zero angle of 
attack are plotted as a function of the experimental chord- 
height ratio. The uncorrected experimental values c/, 
shown by the crosses, are taken from results reported by 
Fage in reference 17. The theoretically corrected values 
Ctf, indicated by the circled points, were computed from 
equation (67) for M'=0. The extrapolated free-air value 
given in reference 17 for each of the bodies is indicated by 
a horizontal dashed line. It is seen that the corrected 
points are in good agreement with the extrapolated free-air 
values. In view of the assumptions made in the theoretical 
development, the relative accuracy of the corrections at 
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lai^e chord-height ratios and lai^e drag coefficients is re- 
markable, particularly in the case of the circular cylinder. 
Glauert (reference 2, pp. 56-57) suggests for the drag 



correction in an incompressible fluid a formula Tvhlch may bo 
written in the notation of this paper as 



l-2Ao— 2, 



<-D(l)l 



(91) 



where (i/c) is the thickness ratio of the airfoil. In this 
equation, as in equation (67), the first correction term appears 
as a result of the interference bet\veen the airfoil thickness 
and the tunnel walls and is identical with the coiTesponding 
term in equation (67) for ikf'=0. The remaining term is an 
empirical correction for the effect of the wake. Tho empiri- 
cal factor K is given by Glauert as a function of (c/0, the 
values being derived by fitting equation (91) to tho experi- 
mental data of reference 17. This wake term differs funda- 
mentally from the wake correction of equation (67) in that 
the correction in this case consists of a single term which 
varies as (c/A) ; whereas the correction in equation (67) com- 
prises two terras, one of which varies as (c/A) and one of 
which varies as (c/A)'. Equation (67) gives corrected results 
which agree as closely with the free-ah- values as do the 
results obtained with equation (91). It has the advantage 
that it is generally applicable to all airfoils and docs not 
depend upon the experimental residts of tests of specific 
sections. 

In summary, the corrected data of figures 9 to 13 indicate, 
for the most part, that when the flow is maintained strictly 
two-dimensional, the theoretical corrections for the tmmel- 
wall interference are, for low Mach numbers at least, accurate 
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FiotntE 13,— Profile drag bg tbree aymmetilcal bodies at lero angle of attaok. 

up to chord-height ratios of unity. The high accuracy 
observed at the larger values of (c/A) must, however, be 
regarded as fortuitous since the theoretical analysis is predi- 
cated upon the assumption that the chord-height ratio is 
small enough that all points of the airfoil may be assumed to 
Uo on the center line of the tunnel and that powers of (c/A) 
higher than the second may be neglected. It is thought 
that, at low Mach numbers, chord-height ratios as high as 
0.7 are permissible if the tests ai'e conducted only for the 
purpose of obtaioing drag characteristics at low values of the 
lift coefficient. However, care must be exercised in ascer- 
taining the maximum chord-height ratio permissible in any 
particular case to insure that the interference lift represented 
by equation (65) is not of such nature and magnitude as to 
affect the general character of the flow in the boundary lay- 
er along the surface of the model. In tests conducted to 
determine the aerodynamic characteristics of a model up to 
and beyond the maximum lift, it is believed that the chord- 
height ratios must be kept to much lower values. At low 
Mach nimibers, chord-height ratios up to 0.4 are probably 
permissible; however, there are no experimental data avail- 
able at present to support this contention. 

As noted previously, no experimental data could be found 
which would permit an evaluation of the accm-acy of the 
calculated effects of compressibility upon the wall-inter- 
ference corrections. Most certainly, as the test Mach 



numbers increase, the permissible chord-height ratios must 
decrease. It is considered that as long as the velocities 
induced at the position of the airfoil by the wall interference 
are small as compared with the velocity of the imdisturbed 
stream, the corrections developed in this paper are appHcable 
even though the stream Mach munber exceeds the critical 
for the airfoil under test. However, as previously noted, at 
and for a range of Mach numbers below choking, the inter- 
ference velocities are no longer small and the corrections are 
iavahd. The extent of this range is imknown. It shoidd be 
emphasized that the flow pattern at and in the immediate 
vicinity of choking does not correspond to any flow pattern 
obtainable with the airfoil in free air; so the test results in 
this range cannot be corrected by any method. 

For zero Mach number (i. e., for an incompressible fluid), 
the results of the present paper can be compared with 
Goldsteiu's particular corrections for airfoils having small 
thickness and camber and small force coefficients. For an 
airfoil on the center line of the tunnel, equations (138), (139), 
(140), (143), and (144) of reference 3, together with the 
expressions of appendix 5 of reference 4, give the following 
equations for the velocity, angle of attack, and aerodynamic 
coefficients in an incompressible fluid: 
F=F'{l+^(2Co-Ci)) 



C,=C,'{1— <r} 

c<,=c/{l-3(7(2(7o-a)} 



(92) 



Here, the moment coefficients are for moments about the 

midchord, and {^^\^ is the moment coefficient at zero lift 

' in free air. The quantities Co and d are determined by 
the shape of the base profile according to the equations 



sm Q 



(7,=- 
" ttJo c si 



cos 2g 
sin B 



dd 



(93) 



In deriving these equations, the notation of references 3 and 
4 has been changed to agree with that of the present paper, 
and the lift-curve slope in free air assumed to have its 
theoretical value of 2ir. 

The corresponding corrections as obtained by setting 
M' = 0 in equations (26), (62), (63), (64a), and (67) of the 
present paper are 

y=y'{l+A<r-f-Tc/} 



cc=a'+, 



c/+4c„/ 



Ci=Ci'{l—<j—2A.ff—2TCa'] 
C<»=c/{ 1— 3Ao— 2tc/} 

;»,=C„/ [l-2Aa-2TC/[ 



(94) 
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The last of these equations is obtained from equations (62) 

' 1 

and (63) by means of the relation c„, =c«. +7 

The correction terms involving a in the two equations for 
the velocity are equivalent, except that the factor A, which 
appears in the equation of the present paper, is replaced in 
Groldstein's equation by the quantity (2(70 — Cj)- Equations 
(93) give 

r2Ci-cy=^ r ^«l:z^d6=« f'^'sin^^ 

IT Jo c sm 9 tJo c 

which becomes, after substitution from equation (44), 

'^«-«=TX'(¥)<f)=i# 

where A is the cross-sectional area of the airfoil. The 
factor A can be expressed in analogous form by means of 

equation (19.05) of reference 2. Since ^ Ac* is equivalent to 

the quantity in reference 2, this equation becomes 



where is the so-called "virtual area" of the base profile. 
The virtual axea of a given body in two-dimensional flow is 
defined as the area occupied by a fictitious quantity of 
fluid having a uniform density p and velocity Vand possessing 
a kinetic energy equal to the total kinetic energy of the field 
of flow about the same body when it is moving forward with 
a steady velocity through an unlimited expanse of incom- 
pressible fluid of density p. The magnitude of the virtual 
area depends upon the shape as well as upon the size of the 
body. It is seen that the first correction term in the velocity 
equation of the present paper (which for the incompressible 
case is simply the result originally derived by Lock) has a 
somewhat higher value than the correction term of the 
Groldstein equation. The Goldstein equation contains no 
term corresponding to the term tc/ in the equation of the 
present paper Groldstein includes this correction, however, 
in the equation for the determination of the stream velocity 
from measurements made at the tunnel wall upstream of 
the model. 

The Goldstein equation for the correction of the measured 
drag coefficient likewise differs from that of the present 
paper by the replacement of the factor A by the quantity 
(2(7o— C^) and by the omission of the term in tc/. If 
Goldstein's equation is apphed to the experimental results 
of Fage given in figiu-es 12 and 13, it is found that there 
is littie to choose between the corrected results given by the 
two equations, except in the case of the circular cylinder 
where the results obtained from the equation of the present 
paper are better. 

The corrections to the lift and moment coefficients as 
derived by Goldstein differ markedly from those of the pres- 
ent paper in that Goldstein's equations contain no terms 



corresponding to the 2A<r and 2rCi' terms which appear in 
the equations of this paper. As has been noted previously, 
the 2tc/ term is accounted for indirectly in the determination 
of the apparent stream velocity. A term of the type 2A(r 
is necessary, however, to correct the measured coefficients 
for the increase in dynamic pressure caused by the inter- 
ference between the walls and the airfoil thickness. 

Since the moment coefficient at zero lift is the same about 
any axis and since the change from the free air to the meas- 
m-ed moment coefficient in the correction to the angle of 
attack win introduce only differences of the second order 
in <r, Goldstein's equation for the corrected angle of attack 
may be written with sufficient accuracy as 

a=«'+2i[./+4(w)J 

In this equation, the part of the correction due to the moment 
on the airfoil is constant, its value depending only upon 
Cm/ for zero lift; whereas in the corresponding equation of 

the present paper the part of the correction due to the mo- 
ment varies with the angle of attack. This difference is of 
small consequence in most applications; however, the equa- 
tion of the present paper, which includes the actual variation 
in moment, may be somewhat the more accurate, especially 
at high angles of attack. 

The compressibility factors which appear in the complete 
equations of the present paper are comparable Avith the 
results of Goldstein and Toxmg (reference 5). The equation 
for drag as given in reference 5, when expressed in coefficient 
form and altered to agree with the notation of the present 
paper, can be written 

c*=c/ { l-[i_(^')»]8/»^<^-[i_(M')»]3/»2«(^) j (97) 

This equation is obtained by modifying equation (91) to 
include the effect of compressibility. Comparison of the 
compressibihty modfficati'ons of equation (97) with those of 
the corresponding terms of equation (67) reveals that the 
compressibility factors appearing in the first correction terms 
differ by the inclusion of a term — (M')' in the numerator in 
equation (67). This difference arises from a failure to note 
in the development of equation (97) that in a compressible 
fluid the dynamic pressure in the timnel is affected by the 
change in density which accompanies the change in axial 
velocity. The compressibility factor of the second (or 
wake-correction) term of equation (97) is not comparable 
with the compressibihty factors of the wake-correction terms 
of equation (67) because of the fundamental difference in 
the nature of the corrections already pointed out in the 
discussion of equation (91). The compressibility factors in 
the equations for lift, moment, and angle of attack in refer- 
ence 5 agree with those appearing in the corresponding 
terms of the equations of the present paper. It should bo 
noted, however, that the lift and moment equations of 
reference 5 include no corrections for the difference between 
the true and apparent dynamic pressures in the tunnel. 
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CONCLUSIONS 

Airfoil data obtained from tests in a two-dimensional-flow 
wind tunnel can be corrected to free-air conditions by means 
of the following equations: 

y=V 1 i+[i-(M')T'^'^+ 1-iM')' J 

2=2 P '[l-(i^')T^ 1-Im'Y 1 ^^^^ 

7?_p/ U 4. 1-0-7 (Mr ■ ^ , [l-0.7(M0'] [l-h0.4(M0^ , 1 

(32) 

[1-1-0.2(MO^[1+0.4(MO^ J 
M=M I 

(33) 



57.3 (T 
'27rVl-(M')* 



(degrees) (64b) 



/fi <r 2-(M0» . 



[2-(M0^[l-|-0.4(M0^] 
1-(M')* 
2-(M0' . 
^~[i-(Ai'0*]*'* 



(62) 



[2-(M0^[l-i-0.4(M0»] , 



+''''4[1-W')*I 



^_ .r 3-0.6 (MQ^ [2- (MQI [1-|-0.4(M01 ^ , 1 



(67) 



where 



and A is a dimensionless factor the value of which depends 
upon the shape of the base profile of the airfoil. (See equa- 
tion (3) and table I.) The remaining symbols are defined in 
appendix B. Nxmierical values of the functions of M' which 
appear in these equations are given in table 11. Eiqperi- 
mental pressure distributions also can be corrected by a 
method outlined in the text. 

The corrections derived are thought to be appHcable up 
to a Mach number near the choking value, which is the 
maximum Mach number attainable in the wind tunnel. The 
choking Mach number is shown to be the stream Mach 
number at which a Mach nmnber of unity is attained locally 
across the tunnel either (1) at the position of the airfoil be- 
cause of the reduction of the available flow area occasioned 
by the presence of the airfoil, or (2) downstream of the airfoil 
as a result of the influence of the airfoil drag upon the flow 
in the wake. The choking Mach number can be estimated 
by means of equations presented in the report. 



Insofar as can be ascertained from the small amoimt of 
experimental data available, the correction equations are 
applicable at low Mach numbers for values of the chord- 
height ratio (filh) as high as 0.7 if the tests ore conducted for 
the purpose of obtaining drag characteristics at low values 
of the lift coefficient. In tests conducted to determine the 
aerodynamic characteristics of an airfoil up to and beyond 
the maximum lift, it is thought that a chord-height ratio of 
0.4 is permissible at low Mach numbers, although there is no 
experimental evidence to support this contention at present. 
At high Mach munbers the permissible chord-height ratios 
must logically be expected to decrease. In particular, if 
the critical speed is exceeded, it is probable that only very 
small values of (c/A) are permissible. There are at this time 
no experimental data available on this aspect of the problem. 

Comparison of the results of the present paper with those 
of references 3, 4, and 6 reveals certain diflferences as noted 
in the section Discussion. 



Ames AERONAtmcAL Laboratort, 
National Advisory Committee fob Abrona-dtics, 
MoPFETT Field, Calif. 

APPENDIX A 

THE VELOCITY FIELD FOR A SOURCE AND FOR A SYM- 
METRICAL AIRFOIL IN A COMPRESSIBLE STREAM 

The velocity induced at a point in a compressible fluid 
stream by a single fluid source con be found to a first degree 
of approximation by a modification of the method used by 
Glauert (reference 7) for the consideration of a vortex in a 
compressible fluid. To this end, a system of polar coordi- 
nates is introduced. The origiu is located at the source and 
the polar axis extends downstream parallel to the velocity V 
of the undisturbed stream. (See fig. 14.) The resultant 
velocity Z7 at any point L{r,<f) is defined by the velocity 
components w and n parallel and normal, respectively, to the 
radius vector. 




FiouBE 14.— Velocity Induced b7 a aoiirce. 



The condition for irrotational motion requires that at all 
points in the field 

c)(m,) dto 



The equation of continuity is 

where pt is the density of the fluid at any point. 



(Al) 



(A2) 
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The source sixength. (mass flow per unit time) is denoted 
by Q. Then, for any cirde enclosing the source, considera- 
tions of symmetry and continuity, respectively, provide the 
two integral relationships 



and 



(A3) 



(A4) 



The radial and circumferential components of the velodly 
may be expanded in the series 



w=y(cos <l>+p^) 



where A, and B, are functions of <j>. If r is laj^e, it is suffi- 
cient to retain only the first terms of each power series, so 
that 



=y (cos 0-1-4) 



^ J 



n= — Y (sin ^+^^ 



(A5) 



To the first power in (1/r), the square of the resultant velocity 
is 

(^ =2^=1 +1 cos 0+J5 sin 0) 

For reversible adiabatic flow, the local density pi is related to 
the density p of the imdisturbed stream by 



Pi=p| 



1 

r-l 



where M is the Mach number of the undisturbed stream and 
7 is the ratio of the specific heats. Thus, to the first power 
in (1/r), 

p,=pjl_M!(^cos<;.-|-5sin<i.)j (A6) 

The solution is now obtained by inserting values from ex- 
pressions (A5) and (A6) into the fundamental equations. 
Substitutions of (AS) into the equation for irrotational 
motion (Al) requires that A shall be a constant. Substitu- 
tion of (A5) and (A6) into the equation of continiuty then 
gives 

^(1 — ilaP sin* <^) =W{A. cos 2^-1-5 sin 2*) 



which becomes upon integration 



B(l— ikP sin» 0) =-| WA sin 2^-1- C 



(A7) 



where C is a constant. The integral equations (A3) and 
(A4) become, respectively. 



Bd(j>=0 



(A8) 



and 



^£Bmx2<t>(h=2^A(l--^—^ (A9) 



Substitution of the expression for B from equation (A7) 
into the integral equation (A9) gives 

^=2;^[vi-iwp] 

while substitution into equation (A8) shows that 0=0. 
Thus, from equation (A7), 



4-jrpl 



r Afsin2<^ "I 
' L(i -ilf* sin* <t>) VT=it2^J 



The expressions for the velocity components therefore 
become 

-^[vi^] 



w=ycos 0- 



2irpr I 



»=— y sin <t>- 



Q 



2irpr 



t M} sin (j> cos <t> ~| 
(l-ikPsin»<^)Vra4^J 



(AlO) 



For a Mach number of zero these equations reduce to the 
well-known results for a source in an incompressible fluid. 

From equations (AlO) the velocity components u and v, 
parallel and perpendicular, respectively, to the direction of 
the undisturbed stream, are found to be 



Q 



COS (f> 



2irpr LVl— ^(1— iWP sin' 

_ Q f yi— M'sin^l 
^~2Tpr\_l-M'sm*<i>J 



(All) 



The drag experienced by the som-ce can be determined 
by evaluating the integral 

D=—j^{pi COS <l>+piw(wco8 <l)—n sin <t>)}r d<l> 

over any circle enclosing the source. To the accuracy pre- 
viously employed, the pressure at any point is 



Pi=p—^-jr(.A cos <t>+B sin ^^>) 



(A12) 



Insertion of this expression, together with (A5) and (A6), 
into the equation for drag gives finally 



D VQ 



(A13) 



which is the same as for a source in an incompressible fluid. 
It is apparent from considerations of symmetry that the Uf t 
force of the source is zero. 

The results of equations (All) can be used to study the 
field of flow about a symmetrical airfoil at zero angle of 
attack in a imiform stream. Such an airfofl can be repre- 
sented by a suitable system of sources and sinks distributed 
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continuously along the chord line. If the notation of figure 
15 is used, the vertical velocity Vi induced in an incompres- 
sible stream, at a given point {xo, yo) on the surface of the 
airfoil is 




FlODBE lEw— Velocity Indnced by a symmeblcal alrfoD. 

whore (j^^ strength per unit length of the source- 

sink distribution in an incompressible stream. From the 
second of equations (All) it follows that the velocity at 
the same point in a compressible stream is 

*'«-2^Jo r(l-Af*sin*,^)^^"S"y 



or 



_2iiEM r 



2kp 



where (^^^ is the strength of the source-sink distribution 

in the compressible case. For any given airfoil of small 
thickness the condition that the flow shall be tangential to 
the surface of the airfoil requires that Vc=Vi at aH points on 
the surface. This fact can be used to relate the source-sink 
distributions for a thin airfoil in the compressible and incom- 
pressible streams by considering the limiting forms of equa- 
tions (Al4) and (Al5) as approaches zero. 

Consider first the limiting form of equation (Al4), which 
may be written 



* 2fl-p Jo (iKb— a;)*+yo* \dx J 



(A16) 



It is seen that even for i/o=0 the integral in this equation is 
finite when evaluated over any interval of integration not 
including the point x=a^. In the limit, the contribution of 
such intervals to the right-hand side of the equation is 
therefore zero, and the equation may be written 



t>t=limlim ^ f^' — 



In evaluating the limit in this equation, care must be taken 
that the limit with respect to 2/0 is taken first in every case. 
Integration by parts gives 



»,=limlim^ 



tan 



«>+« 



+ 



yo 



(A17) 



By virtue of the first mean value theorem for integrals (ref- 
erence 18, p. 65) the integral term iu this equation may be 
written 



tan 



-1 Xq—X 



dx 



where (xo— e) and ^"^^^(.Xo+e). The division into 
two integrals is necessary to ensure that the conditions under 
which the mean-value theorem is applicable are fulfilled: 

namely, that (^^^ is continuous and that tan"' l^as 

the same sign throughout the interval of integration. 
Integration gives 

OT)--^-=[(m..-(m.Jx 
D'--~|io.(>+,^)] 

In the hmit, the value of the terms in the second bracket is 
zero. Thus, only the first term need be retained in equation 
(Al7), which may now be written 

.,=limlimJ-tan-^r(^) +m 1 
= limHmtan-»-i 

Trp\aX /x-xo »->0 iro-HI yo 

Thus the limiting form of equation (Al4) becomes finally 

where and (~^^ ™ay now be considered as pertaining to 

the same general chordwise station z. 

The limiting form of equation (Al5) can be found in 
similar fashion. In this case integration by parts gives in 
place of equation (Al7) 



.=limlimiiziM! 



e-M) Jt-*0 2tP 



io+« 



-1- 



tan" 



dx 



(A19) 



As before, the value of the integral term in this equation is 
zero. The limiting form of equation (A15) becomes finally 



'2p V dx . 



(A20) 



which is the same as (AI8). 
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Since for any given airfoil at all chordwise stations, 

it follows from (Al8) and (A20) that 



(A21) 



that is, the soxirce-siTik distributions necessary to represent 
any given thin symmetrical airfoil in. a uniform stream are 
identical for the compressible and incompressible case. 

This resiilt can be used to calculate the effect of compres- 
sibility upon the field of induced velocities at a considerable 
distance from the airfoil. The increase in longitudinal 
velocity at a large distance yi directly above or below the 
midchord point of the airfoil in an incompressible fluid is 
approximately 

By virtue of the first of equations (All), the corresponding 
velocity at the same point in a compressible fliud is approxi- 
mately. 

Thus, in view of equation (A21), 



(A22) 



that is, in a compressible fluid the increase in longitudinal 
velocity at a point a considerable distance directly above or 
below a symmetrical airfoil is 1/(1— iW*)*" times liie increase 
in longitudinal velocity at the same point in an incompress- 
ible fluid. 

The foregoing results can be used also to determine the 
effect of compressibihty upon the drag of an airfoil in a stream 
having a longitudinal pressure gradient. Consider an un- 
disturbed nonuniform stream having at some given point 
a velocity V, a density p, and a streamwise pressure gradient 
dp/dx. By virtue of Beomoulli's equation, there must be at 
this point a velocity gradient 



dV_ 



(A23) 



This holds true both in the compressible and the incompress- 
ible case. The velocity v a small distance x from the point 
in question is then 



(A24) 



As a result of equation (Al3), the drag experienced by an 
airfoil placed at this point in the stream is, for both the com- 
pressible and incompressible cases, 

where (dQjdx) is, as before, the strength per unit length 
of the somce-sink system necessary to represent the airfoil. 



In order to fulfill the condition that the airfoil is a closed 
body, the soxu:ce-sinfc system must be such that 



Thus the drag is finally 



1 dp 

W2x 



dx 



(A26) 



in both the compressible and the incompressible case. If 
the streamwise pressure gradient is small, equation (A21) 
is still applicable; that is, the source-sink distributions 
necessary to represent the airfoil in the compressible and 
incompressible cases are identical. It therefore follows 
from equation (A25) that the drag of an airfoil in a stream 
having a longitudinal pressiure gradient is unaffected by 
fluid compressibility. 

APPENDIX B 

LIST OP IMPORTANT SYMBOLS 



e 
t 
h 
A 



a 
Ci 

3 

Cd 

V 
M 
M', 
B 



H 
P 
2 
p 

T 
V, 

l+v 

D 

AD 

P 

P. 



airfoil chord, 
airfoil thickness. 

tunnel height. 

a factor depending upon shape of base profile. (See 
equation (3) and table I.) 

'> factor depending upon size of airfoil relative 
to timnel. 

i^^^; factor depending upon size of airfoil relative to 

tunnel, 
angle of attack, 
section lift coefficient, 
section quarter-chord-moment coefficient. 

section midehord-moment coefficient. 

section drag coefficient. 

stream velocity, 
mach number. 

apparent Mach nimiber at choking. 
Keynolds nimaber. 

■^'t ratio of specffic heat of gas at constant pressure 

to specific heat at constant volume (for air 7=1.4). 
total head, 
static pressure. 
djTiamic pressure, 
mass density, 
coefficient of viscosity, 
absolute temperature, 
speed of sound. 

compressibifity factor. (See equation (71) and fig. 4.) 
section drag. 

section drag due to streamwise pressure gradient, 
chordwise lift distribution in coefficient form, 
interference lift distribution. (See equation (78) and 
table ni.) 
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Pj local pressure coefficient. (See equation (69).) 
Si local pressure coefficient. (See equation (68) .) 
X coordinate of points on chord line as measured from 
leading edge. 

6 angxilar coordinate of points on cihord line. (See equa- 
tion (43).) 

r radial distance in polar coordinates. 
(f> polar angle in polar coordinates (positive counter- 
clockwise). 
tp projected thickness of airfoil. 
te effective thickness of airfoil. 
l/i ordinate of base profile. 

^ slope of mean-camber line. 

Q source strength. 

dr 

"Ex ^^^^^'^^'y '^^^ Isiigfcli of chord line. 

u horizontal component of velocity. 
V vertical component of velocity. 

n circumferential component of velocity in polar coor- 
dinates (positive counterclockwise). 

w radial component of velocity in polar coordinates. 

A geometrical area of airfoil section. 

A, virtual area of airfoil section. 

A' cross-sectional area of empty tunnel. 

Am minimum cross-sectional area between model and 
tunnel walls. 

Ai local cross-sectional area of stream tube. 

Ak Fourier coefficients. (See equations (45) and (49).) 

Superscripts: 

(') when pertaining to fluid properties, denotes values 
existing in tunnel far upstream from model; when 
pertaining to airfoil characteristics, denotes values in 
tunnel, coefficients being referred to apparent dy- 
namic pressure g'. 

(") denotes fluid properties far downstream from model. 

(•) denotes airfoil characteristics in tunnel as coefficients 
referred to true dynamic pressure g. 

Subscripts: 

c denotes values ia compressible fluid (excepting V^. 
i denotes values in incompressible fluid. ' 
I denotes local conditions at point in fluid. 
8 denotes conditions existing far downstream when air- 
foil and wake are replaced by source. 



m denotes conditions at miniTmiTn cross-sectional area 

between airfoil and tunnel walls. 
L denotes values on lower surface of airfoil. 
U denotes values on upper surface of airfoil. 
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TABLE I.— VALUES OF A FOR VARIOUS BASE PROFILES 



i— O.XX 


HanVlne 
oval 


Ellipse 


JoukowsU 
section 


Conven- 
tional 
NACA 

sections— 
OOXX 


NACA low-disg sections 


07-oxx: 


lO-OXX 


18-OXX: 


19-OXX 


27-OXX: 


36-oxx: 


46-0X3C 


63-OXX 


64-OXX 


66-OXX 


66-oxx: 


67-OXX 


'o.oo 

.09 
.12 
.IS 
.18 
.21 
.25 
.30 
.35 
.50 
1.00 




0.127 
.106 
.269 
.345 
.426 
.608 
.625 
.780 
.945 
1.500 
4.000 




am 

.173 
.237 
.306 
.376 
.450 
.554 


0.125 
.190 
.284 
.342 
.426 
.512 
.632 


a 119 
.183 
.253 
.326 

.m 

.484 
.696 


ai»4 

.189 
.283 
.341 
.424 
.511 
.631 


ai27 

.196 
.271 
.353 
.439 
.629 
.061 


0.U9 
.186 
.256 
.330 
.408 
.489 
.602 


0.107 
.169 
.217 
.281 
.348 
.417 
.612 


0.106 
.157 
.214 
.278 
.344 
.416 
.513 


a 101 
.148 
.201 
.267 
.316 
.378 
.460 


a 106 
.158 
.215 
.276 
.338 
.402 
.487 


a 109 
.163 
.221 
.283 
.348 
.414 
.601 


a 117 
.174 
.236 
.301 
.368 
.438 
.628 


0.119 
.178 
.243 
.311 
.379 
.450 
.6t3 


a236 
.320 
.403 
.493 
.580 
.703 
.804 
1.019 
1.690 
4.000 


a 165 
.212 
.273 
.337 
.404 
.497 
.626 
.767 

1.258 













































































































BBPORT NO. 782 NATIONAL ADVISORT COMMITTEE FOH AERONAUTICS 



Table H— COMPRESSIBILITY FACTORS FOR CORRECTION EQUATIONS 
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1.029 


1.067 


2.074 


.300 


1.048 


1.0 
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1.152 


1.079 


L173 


1.194 


2.200 


3.394 


1.067 


1.160 


2.174 


.400 


1.091 


1.191 


L2e9 


1.164 


1.341 


L382 


2.390 


3.772 
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1.307 


2.331 


.GOO 


1.155 


1.333 


L&40 


1.270 


1.617 


1. 


m 


2.694 


4. 


!88 


1.210 


1. MO 


2600 


.560 


1.197 


1.434 


L717 


1.353 


1.821 


1.925 


2.914 


4. 


S39 


1.287 


1.706 


2.728 


.600 


1.260 
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1.963 
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2091 
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1.337 


1.910 


2.931 


.625 
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1.627 
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1.379 
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3.0M 
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2279 
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2.4n 
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6,263 


1.428 


2.196 


3. 193 


.673 


1.366 


1.837 


2. 


189 


1.6 


96 


2.716 


2.943 


3.846 


6.788 


1.479 


2.: 


169 


3.364 


.700 


1.400 


1.961 


2 748 


1.804 


3.015 


3.284 


4.148 


7.430 


1.641 


2 676 


3.641 


.725 


1.463 


2.103 


3.061 


1.936 


3.382 


3.701 


4.613 


8.217 


1.613 


2.819 


3.701 


.760 


1.612 


2.2SS 


3.466 


Z095 


3.846 


4.234 


4.1 


168 


9.201 


1.607 


3.115 


4.026 


.776 


1.682 


2.501 


3.982 
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4.433 
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.910 
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97 
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